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Preface 


The  purpose  of  this  study  was  to  compare  several 
techniques  for  the  prediction  of  mass-loaded  natural 
frequencies  and  mode  shapes.  Of  special  interest  in 
this  study  was  the  recovery  of  the  unloaded  mass,  stiff¬ 
ness,  and  damping  matrices  from  measured  modal  data 
using  a  non-square  modal  matrix  and  subsequent  solution 
for  mass- loaded  modal  data. 

This  study  was  somewhat ■ limited  in  scope  in  that 
only  one  panel  with  three  discrete  mass  loadings  was 
experimentally  tested  and  analysed.  It  is  hoped  the 
technique  using  pseudoinverses  will  be  explored  further 
to  determine  the  general  validity  of  this  method.  Appen¬ 
dices  A,  B,  and  C  should  be  helpful  in  this  endeavour. 

I  would  like  to  thank  my  advisor,  Capt.  H.  C,  Briggs 
of  the  Air  Force  Institute  of  Technology,  Mr.  R.  D.  Talmadge 
of  the  Air  Force  Flight  Dynamics  Laboratory,  and  Dr. 

P.  W.  Whaley  of  the  University  of  Nebraska  for  their 
support  and  guidance  in  this  effort.  Additionally,  I 
would  like  to  thank  my  wife,  Janice,  for  her  constant  en¬ 
couragement  and  inspiration.  Finally,  and  most  of  all, 

I  would  like  to  thank  my  Savior,  Jesus  Christ,  i or  that 
"Peace  which  passeth  all  understanding  (Phillipians  -4:7)." 

In  keeping  with  Proverbs  3,  versus  5  and  6, 
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Abstract 


The  purpose  of  this  investigation  was  to  compare  the 
results  obtained  from  three  modal  prediction  techniques . 
The  first  technique  was  an  algorithm  developed  by  Whaley 
for  lightly  damped  structures  (Method  1).  Results  using 
this  algorithm  were  extracted  from  a  thesis  by  Glenesk. 

The  second  method  was  the  finite  element  method  using 
NASTRAN  (Method  2).  The  final  method  was  the  recovery  of 
unloaded  mass  and  stiffness  matrices  from  the  general  ma¬ 
trix-vector  differential  equation  of  modal  analysis  using 
modal  data  obtained  from  an  unloaded  test  item  (Method  3). 
Once  these  matrices  had  been  recovered,  a  quantity  of  mass 
was  added  to  the  mass  matrix  to  simulate  a  mass-loaded 
case.  The  generalized  eigenvalue  problem  was  solved  for 
mass-loaded  frequencies  and  mode  shapes  which  were  com¬ 
pared  to  experimental  results  for  the  same  test  item. 

Both  square  and  rectangular  modal  matrices  were  consider¬ 
ed  in  Method  3.  The  same  test  item  and  three  discrete 
mass-loaded  configurations  which  Glenesk  used  were  tested. 
Percentage  frequency  deviations  from  the  unloaded  test 
item  to  the  mass-loaded  predictions  ranged  from  -1.2-% 
to  +7,4-%  in  Method  1,  from  -20.3-?6  to  +17,84-%  in  Method 
2,  and  from  -20.6-%  to  +8.4-%  in  Method  3.  Several  dis¬ 
crepancies  in  each  technique  prevent  a  direct  comparison 
of  these  results.  The  most  noteworthy  discrepancy  was 


xi 


the  fact  that  the  modal  measurement  procedure  generated 
nonorthogonal  modes.  The  first  method  assumed  the 
mode  shapes  to  be  unaltered  between  the  unloaded  and  mass- 
loaded  cases  while  the  second  method  generated  mutually 
orthogonal  modes.  The  unloaded  nonorthogonal  mode  vectors 
were  used  in  Method  3  to  generated  mass-loaded  modal  quan¬ 
tities.  Detailed  procedures,  results,  and  conclusions 
are  obtained  in  the  body  and  appendices  of  the  report. 
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A  COMPARISON  OF  VARIOUS  TECHNIQUE}  FOR 
THE  PREDICTION  OF  MASS -LOADED  MODE 
SHAPES  AND  NATURAL  FREQUENCIES 


I  Introduction 


Background 

The  ever-expanding  performance  envelopes  of  today's 
highly  complex  fighter  aircraft  subjects  them  to  increas¬ 
ingly  severe  vibration  environments.  Coupled  with  these 
severe  vibration  environments  is  a  desire  to  rapidly  in¬ 
corporate  newly-developed  weapons  system  technologies  into 
the  existing  fleet  of  fighter  aircraft.  One  such  example 
is  the  application  of  laser  physics  technology  to  vibration- 
sensitive  electro-optical  equipment  which  would  subsequent¬ 
ly  be  installed  in  high-performance  fighter  aircraft. 

The  installation  of  electro-optical  equipment  in  an 
aircraft  presents  a  complex  design  problem  in  that  it  is 
necessary  to  know  the  post-installation  modes  of  vibration 
and  natural  frequencies  of  the  aircraft  prior  to  the  actual 
installation  of  this  hardware.  Frequently,  the  only  modal 
data  available  to  the  designer  are  the  pre-installation 
modal  data.  Thus,  due  to  the  vibration  sensitivity  of 
the  electro-optical  hardware,  the  designer  must  consider 
how  to  properly  utilize  the  pre-installation  vibration 
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data  to  correctly  predict  the  post- in.  !  Elation  modes  o: 
vibration,  damping,  and  natural  frequencies.  Whaley  (Ref 
12)  summarized  three  analytical  techniques  to  accomplish 
this  task.  Additionally,  Glenesk  (Ref  3)  utilized  an 
algorithm  developed  by  Whaley  for  lightly  damped  struc¬ 
tures  (Ref  13)  to  predict  the  influence  of  added  lumped 
masses  on  the  vibration  characteristics  of  unloaded  struc¬ 
tures  . 

According  to  the  research  conducted  by  Glenesk  (Ref 
3),  as  the  size  of  the  added  mass  increased  relative  to 
the  mass  of  the  unloaded  structure,  the  accuracy  of  Whaley's 
algorithm  was  signif icantly  affected.  This  degradation 
in  algorithm  performance  might  be  attributed  to  the  assump¬ 
tion  that  the  unloaded  mode  shapes  are  unaffected  by  the 
addition  of  the  lumped  mass.  Another  possibility  involves 
the  fact  that  the  effect  of  damping  was  ignored  in  this 
algorithm.  Thus,  an  added  mass  might  have  significantly 
contributed  to  the  overall  structural  characteristics  in 
such  a  way  as  to  modify  the  mode  shapes  and  natural  fre¬ 
quencies  . 

In  addition  to  the  various  analytical  techniques,  a 
numerical  technique,  finite  element  analysis,  has  been 
widely  used  to  predict  mass-loaded  natural  frequencies, 
damping  ratios,  and  mode  shapes.  This  technique  requires 
the  construction  of  a  computer  model  in  which  the  continuous 


structure  is  idealized  as  a  combi natio'  a  finite 


number  of  various  structural  components  (i.e.,  beams, 
rods,  plates,  etc.).  Although  accurate  results  can  be 
obtained  using  the  finite  element  method,  one  problem 
with  this  technique  is  the  significant  expenditure  of 
human  and  computer  resources  necessary  to  build,  debug, 
and  run  the  finite  element  code. 

The  rapid  development  of  portable  modal  analysis 
equipment  has  made  it  possible  to  eliminate  the  con¬ 
struction  of  a  finite  element  model  entirely.  With 
this  equipment  one  can  lay  out  a  suitable  grid  on  the 
portion  of  structure  to  be  modified,  conduct  standard 
modal  analysis  tests,  and  reduce  the  data  so  obtained 
to  determine  the  desired  unloaded  modal  data.  The 
question  then  becomes  how  to  properly  use  this  data 
to  determine  the  mass-loaded  modal  quantities  for  the 
modified  structure.  One  approach  to  this  dilemma  has 
been  suggested  by  Briggs  and  Whaley  (Ref  1)  whereby  one 
uses  the  general  matrix-vector  differential  equation  of 
structural  analysis,  its  solution  using  a  generalized 
coordinates  approach,  and  the  resulting  definitions 


[u]t[m]  [u] 
[u]t[k]  [u] 

[U]T[C]  [u] 


=  [I] 


2 

(r'n 


n  n 


(1) 

(2) 

(3) 


to  analytically  determine  the  mass-lc  .  >1  modal  data 

when  only  the  experimentally  determined  unloaded  modal 
data  are  known. 

At  the  present  time,  the  usual  solution  to  equa¬ 
tions  1  through  3  requires  that  the  modal  matrices,  [U] 
T 

and  [U]  ,  be  square  matrices.  If  the  mass,  damping, 
and  stiffness  matrices,  [m],  [C],  and  [K],  are  n  x  n 
matrices,  [u]  and  [u]  must  also  be  n  x  n  matrices 
(where  "n"  is  the  number  of  measurement  points).  If 
one  measures  fewer  than  "n"  modes  in  the  frequency 
range  of  interest,  he  must  either  extend  this  frequen¬ 
cy  range  to  accomodate  "n”  modes,  or  reduce  the  grid 
size  to  "n"  grid  points.  As  either  of  these  approaches 
may  be  undesirable,  a  third  approach  using  the  method 
of  pseudoinverses  introduced  by  Penrose  (Ref  10)  may  be 
used  to  isolate  the  mass,  stiffness,  and  damping  ma¬ 
trices  on  the  left  hand  side  of  equations  1  through  3, 
respectively.  Appendix  A  contains  a  sample  problem 
for  the  reader  who  is  unfamiliar  with  this  technique. 
The  resulting  solution  will  be  an  approximate  solution 
to  the  mass,  damping,  and  stiffness  matrices  for  the 
unloaded  structure.  Then,  to  find  the  corresponding 
matrices  for  the  mass-loaded  structure  one  would  add 
appropriate  mass,  damping,  and  stiffness  quantities  at 
the  proper  locations  in  their  respective  matrices  to 
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simulate  the  structural  modification,  .u:d  resolve  equa¬ 
tions  1  through  3  for  the  mass-loaded  modal  information. 
If  damping  is  not  a  factor  one  wishes  to  consider  in 
this  analysis,  one  need  only  consider  the  solution  to 
the  standard  eigenvalue  problem 

[K]  -  oj2  [M]  =  [0]  (4) 

to  determine  the  mass-loaded  natural  frequencies  and  mode 
shapes . 

Purpose 

The  purpose  of  this  investigation  is  to  obtain  data, 
and  compare  the  modal  data  obtained,  using  several  modal 
prediction  techniques.  The  techniques  chosen  for  this 
comparison  are  those  discussed  previously,  namely:  (1) 
GlenesK's  use  of  Whaley’s  algorithm  (Ref  3);  (2)  the 
finite  element  method;  and  (3)  the  method  suggested  by 
Briggs  and  Whaley  (Ref  1).  Since  a  basis  for  comparison 
is  needed,  the  results  obtained  from  modal  prediction 
software  developed  by  Brown  (Ref  2)  will  be  used  as  a 
datum  in  error  precentage  calculations  for  methods  2 
and  3  presented  later  in  this  report.  However,  it  is 
felt  by  the  author  that  to  recalculate  error  values  for 
Method  1  which  would  be  based  on  a  different  datum  would 
be  unfair  to  both  Glenesk  and  Whaley's  algorithm.  Thus, 
all  values  presented  in  reference  to  the  use  of  Whaley's 
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algorithm  will  be  directly  extracted  from  Reference  3. 
The  same  complex,  rib-stiffened  panel  and  several  of 
the  discrete  mass-loading  configurations  investigated  by 
Glenesk  (Ref  3)  will  be  used  in  this  comparison. 

Ob  iectives 

The  objectives  of  this  investigation  are: 

(1)  Experimentally  measure  the  natural  fre¬ 
quencies  and  mode  shapes  of  a  complex  test  panel  in 
the  frequency  range  from  0-500  hz . 

(2)  Construct  finite  element  models  of  this 
panel  and  conduct  a  modal  analysis  on  both  unloaded  and 
mass-loaded  configurations. 

(3)  Use  equations  1,  2,  and  4  along  with 
unloaded  experimental  data  to  obtain  mass-loaded  nat¬ 
ural  frequencies  and  mode  shapes  for  the  test  panel. 

(4)  Present  a  comparison  of  the  results  of 


Objectives  1  through  3 . 


II  Modal  Analysis  and  T q s t  Pro>  •  i  ~os 

Tast  Item 

The  test  item  was  a  panel  fabricated  using  draw¬ 
ings  of  an  upper  fuselage  panel  of  a  C-140  aircraft 
(Figure  1).  The  curved  panel  consisted  of  the  follow¬ 
ing  components : 

1 .  An  outer  skin 

2.  Five  longerons  of  two  different  cross  sections 

3.  Two  curved  main  frame  ring  segments 

4.  Four  edge  doublers,  and 

5.  Various  attachment  hardware  and  bonding  to 
maintain  structural  integrity. 

Mass-Loading  Configurations 

Based  on  the  worst  case  errors  presented  by  Glenesk 
(Ref  3)  for  the  mass-loaded  panel,  four  test  configur¬ 
ations  were  chosen  for  comparison  of  the  three  methods. 
These  were  the  unloaded  panel  and  Glenesk' s  mass- 
loaded  configurations  2 ,  6,  and  7  (Table  1  and  Figure  2). 
The  unloaded  panel  was  included  as  a  means  of  compar¬ 
ing  the  change  in  mode  shape  with  natural  frequency 
which  occurred  between  the  unloaded  case  and  each 
mass-loaded  configuration.  Plots  of  unloaded  versus 
mass-loaded  mode  shapes  allowed  visualization  of  this 
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Panel,  USAF  Drawing  X704D33. 


Table  1 


i  o,  i d  )  MCI 


■  ■ 

Configuration 

Mass  (lb) 

xo  <ft) 

fgna 

2 

0.2420 

1.168 

1 .  168 

6 

0.4158 

1.667 

1.749 

7 

0.1144 

0.499 

0.915 

Figure  2.  Mass- loading  locations  (after  Glenesk 


rx 


0.279 


BAY  1 


■  -  BAY  4 


Notes : 

1. 

Dimensions  in  foot. 

2. 

Jy  specifies  configuration  7 
tion. 

mass- loading  ioca- 

3. 

BAY  1  specifies  internal  un 

supported  bay. 

chanqo  in  mode  shape  ( Figure'  i)  . 

Structural  Models 

Three  grid  sets  were  chosen  to  model  this  panel . 

The  first  grid  set  was  utilized  in  an  effort  to  verify 
Glenesk's  results  by  using  Glcnesk's  grid  set  and 
modal  analysis  procedure.  Glenesk's  model  did  not 
consider  the  discrete  components  as  separate  members. 
Instead,  it  accounted  for  the  total  panel  mass  and 
smeared  this  mass  over  the  grid  which  was  inset  some¬ 
what  from  the  panel  edges.  The  result  was  a  homogene¬ 
ous,  constant  thickness  flat  plate  with  25  grid  points. 
This  model  will  be  referred  to  as  the  Smeared  Stiffener 
Uniform  Model  (Figure  4 ) . 

The  second  grid  was  chosen  to  coincide  with  a 
finite  element  model  which  accurately  modelled  the 
discrete  structural  components  by  allowing  for  panel 
curvature,  discrete  member  cross-sectional  geometry, 
and  offsets  of  component  neutral  axes.  This  model  also 
consisted  of  25  grid  points.  It  will  be  referred  to 
as  the  Discrete  Stiffener  Model  (Figure  5). 

The  final  model  was  chosen  to  investigate  the 
modes  of  vibration  of  each  of  the  internal  bays  (Figure 
2  ) .  This  grid  was  necessary  because  the  previous 
two  models  basically  ignored  the  motion  of  the  unsupported 
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internal  bays.  This  model,  which  con  ;  i  :ed  of  65 
grid  points,  will  be  referred  to  as  the  Bay  Modes 
Model  (Figure  6).  This  model  contained  the  previ¬ 
ous  two  models  as  subsets. 


Modal  Analysis  Test  Procedures 


The  Smeared  Stiffener  Uniform  Model  was  tested 


first.  The  Modal  Assurance  Criterion  (MAC)  developed 
by  Brown  (Ref  3)  along  with  discrete  Transfer  Function 
data  were  used  as  the  basis  for  identifying  candidate 
frequency  ranges  from  which  the  natural  frequencies 
for  each  configuration  were  determined.  The  MAC 
function  is  defined  as  (Ref  3) 


MAC 


Syr(w)  ^ 

S  £  (  (l'  )  S  yy  (  (l!  ) 


where  Syr  =  the  stable  average  of  the  cross  power 
spectrum  between  two  response  measure¬ 
ment  points 

Srr  =  the  stable  average  of  the  auto  power 
spectrum  of  the  stationary  accelero¬ 
meter  response 

Syy  =  the  stable  average  of  the  auto  power 
spectrum  of  the  moveable  accelero¬ 
meter  response. 

Note  that  the  MAC  Function  differs  from  the  more  commonly 
used  Coherence  Function  in  that  the  two  measurements 
in  question  in  the  former  are  two  responses  to  an  im¬ 
pulse  excitation  whereas  like  quantities  for  the  latter 


ItHH 

mm m 


BE! 


r 


would  be  an  impulse  excitation  i nput  md  the  forced  re¬ 
sponse.  In  both  cases,  the  existence  of  a  mode  is  indi¬ 
cated  by  a  region  of  closely  spaced  frequencies  where  the 
MAC  or  Coherence  Function  is  essentially  equal  to  one. 

Since  the  MAC  and  Transfer  Function  data  were  essen¬ 
tially  identical  with  Glenesk's  results,  these  data  were 
not  reduced  to  obtain  natural  frequencies  and  mode  shapes. 
Instead,  Glenesk's  results  (Ref  3)  will  be  used  in  the 
techniques  comparison.  Table  2  contains  a  summary  of  the 
pertinent  data.  Sample  MAC  and  Transfer  Function  plots 
and  data  are  contained  in  Figures  7,  8,  9,  and  10. 

The  Bay  Modes  Model  was  tested  next  using  modal 
analysis  software  developed  by  Brown  (Ref  2).  This  soft¬ 
ware  allowed  the  user  to  select  any  subset  of  the  model 
being  tested  and  consider  only  the  data  relative  to  that 
subset.  Thus,  it  was  not  necessary  to  repeat  this  test 
for  either  the  Smeared  Stiffener  Uniform  Model  or  the 
Discrete  Stiffener  Model.  Representative  mode  shape 
data  obtained  for  these  models  are  contained  in  Figure 
11. 

Modal  Analysis  Data  Reduction 

To  determine  the  natural  frequencies  of  a  given 
model  and  configuration,  one  grid  point  was  selected 
which  was  believed  to  contain  all  of  the  modes  in  the 
frequency  range  from  0-500  hz.  That  is,  it  was  believed 
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Table  2.  Experimental  Natural  Frequencies  Obtained  Using  Modal  Assurance 
Criterion  (MAC)  Function  (Extracted  from  Ref  3). 


Figure  7.  Sample  Modal  Assurance  Criteria  (MAC)  Plot. 
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that  the  point  did  not  1  i .  •  on  or  no  .  r  •  a  node  lino  for 
any  mode.  Examination  of  the?  real  and  imaginary  parts 
of  the  Transfer  Function  for  that  arid  point  revealed 
the  natural  frequencies  for  that  particular  model  and 
configuration.  Detailed  modal  analysis  and  test  pro¬ 
cedures  are  contained  in  Appendix  B.  Data  acquisi¬ 
tion  programs  for  use  in  the  Hewlett-Packard  HP-ddilB 
Fourier  Analyser  are  presented  in  Appendix  C. 

Once  the  natural  frequencies  for  a  particular  model/ 
conf igura t ion  had  been  identified,  the  Transfer  Function 
data  from  each  point  on  the  grid  in  question  were  reduced 
to  yield  mode  shape  vectors  for  each  mode.  For  purposes 
of  comparison  with  the  finite  element  method  the  discrete 
mode  shape  vectors  were  converted  to  a  format  identi¬ 
cal  to  NAS TRAN  output.  NASTRAN  data  wore  processed  usinq 
a  standard  graphics  package,  HCSXAST  (Ref  6) ,  which  was 
used  to  display  plots  of  the  undeformed  versus  deformed 
mode  shapes  in  both  unloaded  and  mass-loaded  conf iqurat lot 

Resul ts 

Modal  Assurance  Criteria  Function  test  results  are 
presented  in  Table  2.  Corresponding  results  from  the 
Modal  Analysis  testing  are  shown  in  Table  3 .  Compar isan 
of  the  data  in  these  two  tables  reveals  numerous  areas 
of  disagreement  in  the  natural  frequencies  of  each  con¬ 
figuration.  This  apparent  discrepancy  can  be  explained 
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when  one  considers  the  tost  toch.n i gu'* 
set  of  data  were  obtained.  The  MAC  testing  used  response 
data  from  a  fixed  reference  accelerometer  and  a  moveable 
accelerometer,  and  15  impulse  excitations  which  were 
randomly  spaced  over  the  entire  panel.  The  spatial 
randomness  of  the  excitations  enhanced  the  probability 
that  every  mode  in  the  structure  would  be  excited  in 
that  not  every  excitation  would  lie  on  a  node  line.  In 
contrast,  the  Modal  Analysis  testing  used  a  fixed  ex¬ 
citation  point  in  conjunction  with  a  moveable  accelero¬ 
meter  to  measure  the  structural  forced  response.  Inherent 
in  the  latter  technique  is  the  assumption  that  the  chosen 
excitation  point  never  lies  on  or  near  to  a  node  line 
and  thus  the  mode  shapes  obtained  from  this  technique 
represent  a  unique  set  of  modal  data.  Therefore,  if 
the  chosen  excitation  point  lies  on  or  near  to  a  node 
line,  some  modes  may  be  "missed"  during  examination  of 
Transfer  Function  data  because  these  modes  were  never 
excited  to  begin  with.  In  this  respect  it  would  appear 
that  the  MAC  function  data  may  bo  the  more  accurate  data. 
Since  a  comparison  of  modal  prediction  techniques  is  the 
thrust  of  this  report  this  comparison  will  be  carried  out 
only  on  those  data  for  which  corresponding  results  in 
the  MAC  Function  data  are  available.  This  comparison 
method  will  be  used  throughout  the  remainder  of  this 
report. 
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Ill  The  W ha  1  o y  Algor  i  f  I  in  Me  * 


Overview 

The  algorithm  developed  by  Whaley  (Ref  13)  and  sub¬ 
sequently  used  by  Glenesk.  (Ref  3)  substitutes  suitable 
expressions  for  the  kinetic  and  potential  energy  of  a 
flat  plate  into  Lagrange's  equations  of  motion.  Then, 
taking  the  first  variation  of  the  expression  for  the 
virtual  work  of  the  applied  inertial  loads  with  respect 
to  the  generalized  coordinate,  the  expression  for  the 
generalized  force,  Qj_,  is  obtained.  The  final  form  of 
this  expression  is 
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Equation  5  contains  the  effects  of  the  inertial  forces 


due  to  the  added  lumped  mass  (M0,  Rx,  Rv>  xD,  and  yD) 

which  are  proportional  to  the  second  derivative  of  the 

d2qj 

generalized  coordinate,  - j—  or  q  .  These  effects,  m 

dt^  J 

turn,  become  additions  to  the  mass  matrix  when  the 


eigenvalue  problem  is  solved  for  natural  frequencies 


and  mode  shapes . 
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Result  s 

The  results  from  Whaley's  algorithm  are  summarized 
from  Glenesk  (Ref  3)  in  Tables  4,  3,  and  6.  The  data 
generated  using  the  MAC  Function  wore  merged  with  Whaley's 
algorithm  (Ref  13)  to  yield  the  Whaley  Algorithm  results. 
Glenesk  (Ref  3)  noted  that  a  comparison  of  unloaded  versus 
mass-loaded  mode  shapes  was  the  means  of  determining  the 
unloaded/predicted  mass-loaded  frequency  pairings. 

Examination  of  the  data  in  Tables  4,  5,  and  6  reveals 
the  largest  percentage  error  between  actual  and  predicted 
values  occurred  on  Configuration  7  (Table  6,  Mode  1)  where¬ 
as  the  smallest  percentage  error  value  occurred  on  Con¬ 
figuration  2  (Table  4,  Mode  9).  With  the  exception  of 
Mode  1,  Configuration  7,  all  predicted  frequency  values 
were  within  +  7.5 %  of  the  experimentally  measured  values. 

Table  7  presents  a  comparison  of  the  absolute  per¬ 
centage  in  frequency  shift  from  the  unloaded  panel  to 
each  of  the  three  mass-loaded  conf igurations .  The  largest 
percentage  frequency  shifts  occur  when  the  largest  mass 
was  located  on  an  unsupported  portion  of  the  panel  skin. 
With  the  exception  of  modes  8  and  9  for  Configuration  6, 
all  frequency  shifts  were  within  +  3'Y.  of  the  unloaded  fre¬ 
quency. 
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Table  4.  Data  Results  -  Configuration  <  extracted  from  f 


Testina  Technique 


Percent 


Mode 


1 


Modal  Assurance  Criteria 


( 

174.08 

187.62 

204.00 

225.31 

243.32 
263.80 
282.45 
363.08 


Table  5.  Data  Results  -  Configuration  6  (Extracted  from  Ref 


Mode 


1 


Testina  Technique^ 

Modal  Assurance  Criteria 

Whaley  Alqorithn 

169.71 

166.22 

180.72 

175.34 

199.30 

191.23 

208.27 

205.63 

238.22 

235.20 

252.20 

251.18 

279.95 

278.44 

292.20 

271.14 

360.20 

338.70 

Percen 

error 


-2.06 

-2.98 

-4.05 

-1.27 

-1.27 

-0.40 

-0.54 

-7.21 

-6.11 
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Table  6.  Data  Results  -  Configuration  .  (extracted  from  Ref 

3). 


Mode 

Testing  Techniques 

HEB9SS 

wBg<5>m<5ianaar.BaM«g4»M 

Wha In v  Alaorithm 

'l 

143.10 

167.77 

17.24 

2 

175.09 

180.20 

2.91 

3 

184.83 

191.16 

3.44 

4 

204.30 

206.98 

1.31 

5 

238.22 

237.33 

-0.37 

6 

257.50 

250.45 

-2.73 

279.88 

278.70 

-0.42 

292.76 

292.48 

-0.09 

361.71 

364.39 

0.74 
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Table  7.  Percentage  of  frequency  Shift  Using  Whaley's  Algorithm  versus  Unloaded 
Measured  Data  (Extracted  from  Ref  3). 
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IV  The  Finite  Element  X~ ;  '  -i 


Overview 

Several  excellent  finite  element  codes  for  struc¬ 
tural  analysis  are  presently  in  existence.  One  such 
widely  used  program,  NAS TR Ah'  (NAsa  STRuctural  ANalysis, 
Ref  9,11),  was  selected  for  use  in  the  finite  element 
modelling  and  modal  analysis  of  the  structure.  This 
code  includes  the  general  20  degree-of-f reedom  quadri¬ 
lateral  elements  (CQUAD2)  and  12  degree-of-f reedom  bar 
elements  (CBAR)  of  which  the  test  structure  was  con¬ 
structed.  NASTRAN  also  contains  provisions  to  allow  for 
the  offset  of  the  neutral  axes  of  the  bar  elements  from 
the  grid  points  which  were  defined  at  the  midsurface  of 
the  panel  skin.  Thus,  the  cross-sectional  and  spatial 
properties  of  each  component  stiffener  could  be  included 
in  the  analysis.  Onlv  the  out-of-plane  component  of  the 
vibration  (i.e.  -  radial  component  for  curved  models 
and  z-component  for  the  flat  model)  was  investigated. 

Finite  Element  Models 

Three  finite  models  were  constructed  to  coincide 
with  the  test  grids  described  in  Section  II.  The  Smear¬ 
ed  Stiffener  Uniform  Model  considered  the  panel  to  be 
a  flat  plate  (Ref  3).  This  model  did  not  extend  to  the 
panel  extremities  but  was  inset  somewhat  from  the  panel 
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edges  (Figure  2)  to  coincide  with  the  grid  used  by 
Glenesk  (Ref  3) .  The  mass  of  the  entire  panel  includ¬ 
ing  stiffeners  was  smeared  over  this  grid  to  provide 
a  uniform  thickness  model  with  homogeneous  material 
properties.  Although  the  mass  of  this  model  was  identi¬ 
cal  to  the  overall  structure  mass,  this  model  was  more 
dense  than  had  it  been  extended  to  the  geometric  test 
panel  boundaries.  It  consisted  of  25  grid  points  from 
which  data  were  obtained  and  16  quadrilateral  elements. 

The  Discrete  Stiffener  Model  (Figure  4)  was  con¬ 
structed  to  allow  for  panel  curvature,  discrete  stiff¬ 
ener  geometries  (i.e.  -  different  cross-sections),  and 
the  offset  of  the  stiffener  neutral  axes  from  the  panel 
surface.  This  model  was  designed  to  faithfully  repre¬ 
sent  the  panel  from  a  structural  standpoint  while  main¬ 
taining  the  25  grid  points  of  the  Smeared  Stiffener 
Uniform  Model.  Forty- four  bar  elements,  25  grid  points, 
and  16  quadrilateral  elements  were  used  in  this  model. 

The  Bay  Modes  Model  (Figure  5)  was  included  to  in¬ 
vestigate  the  motion  of  each  of  the  four  internal  bays. 
This  model  was  an  extension  of  the  Discrete  Stiffener 
Model  in  that  two  extra  sets  of  five  grid  points  per 
bay  were  added  to  that  model  to  derive  this  model.  This 
model  was  constructed  of  65  grid  points,  48  quadrilateral 
elements,  and  76  bar  elements. 

Appendix  C  contains  a  sample  of  the  Bulk  Data  Decks 
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which  wore  used  to  qenor.it,'  irod.il  d  c  :  or  t  h«>se  model  . 

Results 

Only  the  Bay  Modes  Model  result.'-;  (Table:-  8  throuqh 
11)  will  be  discussed  since  this  model  yielded  the  be  L 
representation  of  the  overall  panel  motion.  As  in  the 
use  of  the  Whaley  Algorithm,  a  pairing  of  unloaded  versus 
mass-loaded  mode  shapes  (Figure  12)  was  used  to  track  the 
change  in  natural  frequency  between  the  unloaded  panel 
and  each  mass-loaded  configuration.  Upon  observation  of 
Figure  12  it  is  seen  that  this  process  is  somewhat  subjec¬ 
tive  in  nature  in  that  the  mode  shapes  do  not  remain  com¬ 
pletely  unaltered.  It  is  left  to  the  discretion  of  the 
engineer  to  properly  select  the  mode  pairings,  and  herein 
lies  a  potential  source  of  error.  In  an  attempt  to  have 
this  mode  pairing  as  unbiased  as  possible,  the  author  con¬ 
sulted  another  engineer  to  independently  aid  him  in  this 
process.  Mode  pairings  which  were  not  in  agreement  between 
the  author  and  the  other  engineer  were  discussed  and  a 
consensus  of  opinion  arrived  at.  The  results  from  this 
mode  shape  pairing  exercise  are  presented  in  Table  12. 

Here  absolute  changes  in  natural  frequency  from  the  un¬ 
loaded  to  mass-loaded  panel  range  from  essentially  zero-’, 

(Mode  20,  Configuration  2)  to  -20. 3-",  (Mode  16,  Configura¬ 
tion  6).  Out  of  57  such  pairings,  the  frequency  shift 
from  the  unloaded  to  the  mass-loaded  panel  in  was  within 

l 

±  7.5-%  of  the  unloaded  panel  in  70  percent  of  the  pair  inns  j 
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Table  8.  N’ASTKAN  Modal  Data  -  Unloaded  Pane 


Table  9.  .TASTRAN  Modal  Data  -  Conf iqurat i on 


Table  10.  NASTRAN  Modal  Data  -  Configuration  6 


in 


Figure  12.  Mod.'  .il  Pipes,  Bay  Modes  Finite  Kletr.ent  M 
for  A.)  '"lean  Panel  (20G.BG  Hr. )  ,  B.)  Co 
ration  2  (2nB.ll  Hz),  c.)  f’onf  iourai  ion 
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and  within  +  3-°..  of  tho  unloaded  panel  "u  percent  of  th 
time . 


V  The  Use  of  Psoudn  i  nv; ; r  ■  '  •  n  lb'  _  2 )Vl'>i'V  of 

the  Discrete  Mass  ,  S  L  i. f  ness  ,  and  .^a.ioiiiq 
Matrices  and  Sol  nt  i  on  of  L iv ' 

Mass- Loaded  Problem 


Overview 

The  method  suggested  by  Briggs  and  Whaley  (Ref  1) 
solves  the  general  matrix- vector  differential  equation 
using  generalized  coordinates  to  obtain  the  definitions 
of  equations  1  through  3.  Then,  using  modal  data  ob¬ 
tained  from  standard  vibration  test  methods,  and  suit¬ 
able  matrix  manipulations  (See  Appendix  A) ,  equations 
1  through  3  are  solved  for  the  mass,  stiffness,  and 
damping  matrices.  The  usual  method  of  solution  re¬ 
quires  a  square  modal  matrix  which  is  inverted  in  the 
solution  for  these  matrices.  However,  in  general, 
this  matrix  may  be  rectangular  with  more  rows  (grid 
points)  than  columns  (mode  vectors) .  Thus,  one  would 
like  to  have  a  means  to  solve  equations  1  through  3 
when  the  modal  matrix  is  non-square.  The  pseudoinverse 
(Ref  8)  presents  such  a  method  for  inverting  rectangu¬ 
lar  matrices  and  thus  potentially  for  a  more  general 
solution  to  these  equations. 

Data  Reduction 

A  computer  program  was  developed  to  take  the  ex¬ 
perimentally  determined  unloaded  modal  data,  compute 
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the?  unloaded  modal  matrix  and  required  inverses,  and 
solve  for  the  unloaded  mass,  stiffness,  and  damping 
matrices.  Since  only  the  natural  frequencies  and  their 
resulting  mode  shapes  were  of  interest  in  this  study, 
the  standard  eigenvalue  problem  posed  by  equation  4  was 
solved  with  suitable  additions  of  mass  in  the  mass  matrix 
for  the  mass-loaded  natural  frequencies  and  mode  shapes. 
Appendix  B  contains  a  more  detailed  explanation  of  the 
data  reduction  process;  Appendix  C  contains  a  copy  of 
the  program  used  for  data  reduction. 

Results 

An  extremely  simplified  experimental  set  up  consist¬ 
ing  of  nine  grid  points  was  used  for  each  mass-loaded  con¬ 
figuration  (Figure  13)  to  demonstrate  the  validity  of  the 
computer  program  before  extension  to  the  more  general 
pseudoinverse  case  was  attempted.  The  structure  was  test¬ 
ed  in  both  unloaded  and  three  discrete  mass-loaded  con¬ 
figurations  from  which  nine  frequency/mode  shape  pairs 
were  identified  for  each  unloaded/mass-loaded  configura¬ 
tion  (Figure  14).  The  nine  modal  vectors  were  used  to 
form  a  square  9X9  modal  matrix;  the  nine  natural  fre¬ 
quencies  wore  used  to  form  the  matrix  c n  the  right  side 
of  equation  2.  Then,  using  the  technique  described  in 
Appendix  B,  the  unloaded  mass  and  stiffness  matrices  wore 
recovered,  and  a  quantity  of  mass  equal  to  the  added  mass 
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Figure  14.  Nine  po j. n t  ronrsp  grid  mod'*  snipe-,  eon:  jcjur 
tion  7,  "ode  1,  for  A.)  I'n  loaded  ['.me!  (HVi. 
Hz),  H.)  l.o. id>'d  Panel  (13‘'.')7'.  Hz),  n.)  ;qu,. 

Modal  Matrix  Predict  ion  (  ]  f»H  ,  8-17  Hz)  and  :>,  ) 
Pseudo  inv.-r  :e  Predict.  ion  (  Hi'1.  1  17  Hz).  ( i>»f 
ed  panel  denoted  by  dotted  lines) 
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was  added  to  tho  5,5  term  (i.o.,  tin-  iss  loading  poi.r.;, 
See  Figure  13)  of  the  mass  matrix.  Finally,  equation  4 
was  solved  using  a  generalized  eigenvalue  routine,  ETC-ZF 
(Ref  5),  for  the  mass-loaded  natural  frequencies  and  mod' 
shapes.  These  predicted  mode  shapes  were  then  compared 
to  the  actual  measured  mode  shapes  to  verify  that  tho 
program  had  successfully  predicted  both  the  correct  natu¬ 
ral  frequency  and  its  corresponding  mode  shape.  From 
this  mode  shape  comparison  (Tables  13,  1-1,  and  15)  it 
was  determined  that  the  program  using  a  square  modal  ma¬ 
trix  had  correctly  predicted  the  frequency  within  _+  5-°o 
error  2-1  out  of  27  times,  and  with  ±  10 %  error  26  out 
of  27  times.  These  data  are  presented  in  Tables  13,  14, 
and  15. 

Next,  the  last  column  of  the  modal  matrix  and  tho 
last  row  and  column  of  the  matrix  of  natural  frequencies 
squared  were  deleted  to  simulate  a  case  whore  fewer  than 
"n"  modes  were  measured  where  "n"  is  the  number  of  grid 
points  (or  rows  in  tho  modal  matrix).  The  modified  9  X 
8  modal  matrix  was  then  inverted,  and  9X9  mass  and 
stiffness  matrices  were  generated.  The  mass  matrix  was 
perturbed  by  addition  of  a  quantity  of  mass  equal  to  tho 
mass  loading  conf  i  oration  at  the  5,5  location,  and  tlx.' 
general  eigenvalue  problem  was  again  solved  using  E1CZK. 
Upon  examination  of  the  actual  unloaded  and  predicted 
mass-loaded  psnudninverse  frequency  data. 


it  was 


found 


that  these  data  are  identical  to  six  decimal  places. 

Of  even  more  concern  is  the  fact  that  the  predicted 
mode  shapes  were  quite  different  from  the  measured  mode 
shapes.  Thus,  from  this  preliminary  investigation,  it 
did  not  appear  that  the  method  employing  pseudoinverses 
would  yield  valid  results  and  further  attempts  at  its  use 
were  abandoned. 

However,  it  was  demonstrated  that  the  method  using 
square  modal  matrices  will  yield  valid  results.  One  draw¬ 
back  in  the  use  of  square  modal  matrices  in  this  method 
is  that  when  one  increases  the  dimension  of  the  modal 
matrix  (i.e.  -  the  number  of  grid  points  or  rows),  a 
corresponding  number  of  mode  shape  vectors  and  natural 
frequencies  must  be  generated.  At  the  outset  of  this  re¬ 
port  it  was  pointed  out  that  for  larger  numbers  of  grid 
points  this  may  be  both  undesirable  and  that  it  may  not 
be  possible  to  obtain  a  large  number  of  natural  frequen¬ 
cies.  Thus,  there  is  a  need  to  limit  the  number  of  grid 
points  when  one  models  the  structure  in  question.  Un¬ 
fortunately,  there  is  at  present  no  g onoral  method  to 
predict  the  number  of  natural  f roquoncics/modos  within 
a  given  frequency  ranqe.  The  solution  to  this  problem 
may  be  to  initially  begin  with  a  simplified  nine  point 
grid,  test  the  unloaded  structure  and  reduce'  the'  Trans¬ 
fer  Function  data  to  obtain  the*  numbeu  e>f  natural  f re." 
quencie's/modexs  in  the  specific  frequency  rang*'  of  interest  . 
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This  number  of  natural  f requeue ios/m  k:  will  give  an 

indication  of  the  maximum  number  of  grid  points  avail¬ 
able  if  one  wishes  to  refine  this  grid  and  use  the  square 
modal  matrices.  Since  the  Transfer  Function  data  theoret¬ 
ically  contain  the  same  information  for  any  grid  point 
unless  the  response  accelerometer  for  a  particular  grid 
point  was  situated  on  a  node  line,  it  is  not  required 
to  initially  increase  the  number  of  grid  points  above 
nine.  During  data  reduction  to  determine  the  maximum 
number  of  natural  frequencies  (or  grid  points),  it  will 
thus  be  advantageous  to  reduce  the  data  from  more  than 
one  grid  point  to  insure  that  no  modes  were  "missed"  duo 
to  a  given  accelerometer  being  inadvertently  placed  on 
a  discrete  node  line. 


4R 


VI  Discussion  of  t ho  Motor 


General 

As  with  any  endeavour  of  this  typo,  a  learning  curve* 
is  associated  with  one's  ability  to  use  unfamiliar  equip¬ 
ment,  software,  and  testing/modelling  techniques.  This 
was  especially  true  for  the  author  who  had  no  previous 
experience  in  modal  analysis  testing,  with  the  associated 
data  reduction  techniques,  nor  with  finite  element  model¬ 
ling  procedures.  Thus,  the  author  was  in  a  unique  posi¬ 
tion  to  bo  able  to  evaluate  the  three  modal  analysis  tochni 
ques  which  are  the  subject  of  this  report.  Admittedly 
this  assessment  of  the  three  methods  will  be  only  one 
person’s  viewpoint  and,  as  such,  is  somewhat  subjective 
in  nature.  However,  it  represents  the  viewpoint  of  otic 
who  was  previously  uninitiated  in  this  area  of  expertise. 
The  evaluation  of  each  method  will  bo  presented  separate¬ 
ly  in  the  successive  paragraphs  of  this  section. 

At  the  outset  of  this  investigation  it  was  cited 
that  tiie  Modal  Analysis  Software  developed  by  Brown  (Ref 
2)  would  be  used  as  a  datum  for  comparing  the  methods. 
However,  because  the  mode  shapes  from  Brown's  software 
and  NASTRAN  were  quite  different,  such  a  comparison 
between  experimental  data  and  NASTRAM  data  was  not  possi¬ 
ble.  Furthermore,  since  the  method  using  psoudoima -it  :t 
was  abandoned  in  favor  of  square,  invertible  0  X  o  modal 
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matrices,  and  extension  to  either  2  or  6a-  arid  point..-; 
was  not  possible  due  to  a  lack  of  a  sufficient  number  of 
natural  frequencies  in  the  specific  range  from  0-500  Hz, 
a  comparison  of  this  method  to  Brown' s  software  was  like¬ 
wise  not  possible.  Neither  was  it  possible  to  tie  Glenesk's 
modal  data  to  corresponding  modal  data  generated  by  Brown's 
software  due  to  lack  of  sufficient  modal  data  in  Glenesk’s 
report  (Ref  2).  Thus,  each  method  will  be  discussed  as 
fully  as  possible  while  keeping  these  limitations  in  mind. 

The  Whalev  Algorithm 

The  results  obtained  from  use  of  the  Whaley  Algorithm 
were  extracted  directly  from  Reference  3  and  are  presented 
in  Tables  4  through  6.  Because  the  author  did  not  direct¬ 
ly  use  this  algorithm,  the  author  is  unable  to  present 
an  evaluation  of  the  difficulties  encountered  in  obtain¬ 
ing  results  from  this  method.  However,  when  one  considers 
the  modes  which  were  predicted  by  the  algorithm  when  modal 
testing  did  not  reveal  such  a  mode,  one  can  envision  seri¬ 
ous  limitations  to  its  use.  As  previously  noted,  Olenosk 
reported  a  degradation  in  algorithm  performance  as  the 
size  of  the  added  lumped  mass  increased  relative  to  the 
overall  structure  rasa.  Thus ,  an  as  yet  undetermined 
upper  bound  exist:-:,  beyond  which  the  predicted  mode  :-,n,,pe:, 
and  natural  frequencies  will  cease  to  have  an  accept  able 
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degree  of  accuracy. 


one  to  obtain 


The  use  of  this  a  Igor  i  tin:,  requ 
modal  data  on  tho  unloaded  structure  before  the  mass- 
loaded  modal  quantities  can  be  predicted.  This  requires 
modal  testing  facilities,  test  hardware,  and  the  nccossa 
software  for  conducting  tho  required  Fourier  analysis 
on  the  Transfer  Function  data.  This  can  moan  a  signi¬ 
ficant  commitment  in  terms  of  resources  (manpower,  money 
testing  facility,  etc).  However,  with  the  advent  of 
portable  modal  analysis  equipment,  this  investment  in 
resources  is  expected  to  decline  rapidly.  One  factor 
not  as  easily  evaluated  is  the  level  of  expertise  re¬ 
quired  to  successfully  conduct  the  required  modal  test¬ 
ing  and  subsequent  data  reduction  in  order  to  obtain 
valid  results.  It  was  the  author's  experience  that  one 
must  rely  heavily  upon  "experts"  in  this  field  to  avoid 
the  pitfalls  of  questionable  test  procedures  or  of  on 
invalid  data  analysis.  Tho  knowledge  required  to  suc¬ 
cessfully  conduct  t he  testing  and  data  analysis  is  one 
of  those  intangible  factors  upon  which  it  is  extremely 
difficult  to  place  an  evaluation. 

The  one  serious  limitation  inherent  in  the  use  of 
this  algorithm  is  that  it  will  not  predict  mass-loaded 
mode  shapes .  Inst  >ad,  it  considers  tho  mode  shapes  to 
be  unaltered  by  the  addition  of  mass  to  the  structure. 
Although  unilt'-n  d  mode  shapes  were  observed,  they  were 


randomly  sea tiered  throughout  the  .  :rimont.-.«l  dat  <t . 

It  is  the  opinion  of  the  author  upon  review  of  the  ex¬ 
perimental  data  that  even  with  a  frequency  shift  the 
mode  shapes  were  altered  far  more  often  than  not.  This 
is  not  believed  to  be  a  function  of  experimental  methodo¬ 
logy  or  questionable  data,  but  rather  an  actual  occurrence 
in  nature. 

The  Finite  Element  hot hod 

Again,  as  in  the  case  of  the  modal  analysis  testing, 
a  learning  curve  was  associated  with  the  use  of  the 
finite  element  method  as  presented  in  NT  AS  TRAN .  Unlike 
modal  analysis  testing  where  the  methodology  used  to 
obtain  and  reduce  the  data  is  all  important,  the  author 
found  the  structural  modelling  technique,  especially 
modelling  of  the  boundary  conditions,  to  be  critical  in 
obtaining  good  results.  Since  the  structure  was  support¬ 
ed  by  fcungy  cords,  the  support  conditions  fell  somewhere 
between  the  "free-free"  case  and  the  "clamped-clanpod" 
case  at  the  panel  boundaries.  Thus,  it  was  not  possible 
to  clamp  the  finite  element  model  at  the  odoes ,  nor  was 
it  possible  to  allow  the  panel  to  be  totally  unrest ru ined. 
Several  methods  of  modelling  the  boundary  conditions 
were  used  to  determine  the  model  which  most  closely  ap¬ 
proximated  the  actual  boundary  ci.nJ:  t  ions.  Among  the 

wem  the  restraint  (if  rigid  body  modes  by 
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methods  tried 


modal  l.i.  nq  t  ho  bunqy  cords  Ion;  •  1.  :  !<  r  rods  with 
stiffness,  the  uso  of  the  SUI’ORT  card  or  SPCl  cards  to 
restrain  suitable  deqreus-uf-f roodor,  at  throe  non-col  ino-»r 
qrid  points  on  the  panel  edges,  and  the  use  of  the  SUPOHT 
card  to  restrain  all  six  degree-of-f roedom  at  one  select¬ 
ed  grid  point.  In  the  latter  case  a  COMM  2  card  was  used 
with  small  mass  and  inertia  values  at  the  selected  grid 
point  to  alleviate  the  singularity  in  the  mass  matrix  im¬ 
posed  by  restraining  six  degrees-of-f reedom  at  the  same 
grid  point.  The  technique  using  SPCl  cards  was  chosen  as 
the  best  method  to  model  the  supports. 

In  general,  each  finite  element  model  generated  a 
greater  number  of  natural  frequencies  (Tables  8  through 
12)  than  the  MAC  function  test  results  (Tables  2,  -1,  5, 
and  6)  indicated  wore  present  in  the  frequency  range  from 
0-500  Hz.  However,  when  the  finite  element  results  (Table 
12)  were  compared  to  Modal  Analysis  test  results  (Table 
7),  approximately  the  same  number  of  natural  frequencies 
were  obtained  using  these  two  methods.  Although  the  quani- 
ty  of  frequencies  obtained  using  the  finite  element  method 
was  approximately  the  same  as  those  obtained  from  Modal 
Analysis  results  and  the  frequencies  wore  numerically 
similar,  it  became  apparent  upon  comparison  of  the  mode 
shapes  from  numerically  similar  frequencies  that  more 
numerical  similarity  of  the  natural  f roquencios  was  in¬ 
adequate  by  itself  as  a  criterion  for  comparing  the 
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Modal  Analysis  tost  results  with  sin.il  •  r  results  f ron 
NASTRAN .  As  the  desired  means  of  comparison  of  those  tv;o 
sets  of  data  was  a  "pairing"  of  "like"  mode  shapes,  it 
was  unfortunate  that  this  was  not  possible  due  to  the 
large  discrepancy  in  mode  shape  versus  natural  frequency 
when  the  two  sots  of  data  were  compared.  Thus,  a  better 
means  of  comparing  analytical  to  experimental  data  is 
needed  and  should  be  the  subject  of  further  research. 

This  failure  in  ability  to  compare  experimental  to 
analytical  data  presents  a  perplexing  problem  for  the  de¬ 
signer  who  wishes  to  minimize  the  vibrational  effects  on 
the  hardware  to  be  installed  by  placing  the  associated 
attachment  hardware  at  points  of  minimum  vibration.  Which 
mode  shapes  does  the  designer  believe?  Even  more  distres¬ 
sing  are  the  "additional"  modes  which  were  computed  by 
NASTRAN.  Do  these  modes  really  exist?  Are  they  heavily 
damped  modes  which  exist  not  far  from  the  noise  floor? 

Were  these  modes  missed  in  the  experimental  testing  due 
to  an  inadequate  sampling  bandwidth?  Or  are  these  "addition¬ 
al"  modes  purely  "synthetic"  modes  generated  by  the  solution 
software  in  NASTRAN  but  not  actually  existing  in  nature? 

These  questions  remain  unanswered  and  are  seen  as  a  severe 
limitation  of  the  "state-of-the-art"  in  vibration  analysis. 

Although  a  pairing  of  experimental  versus  analytical 
(finite  element,  NASTRAN-gone rated)  mode  shapes  was  not  suc¬ 
cessful,  no  such  lack  of  success  was  encountered  when 
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pairing  either  XASTRAN-qoneratod  uni  c  e..ed  data  to  XASTRV  - 
generated  mass-loaded  data  (Table  12  and  Figure  12)  or 
when  corresponding  Modal  Analysis  Software  data  (Table  3) 
were  compared.  Thus,  both  methods  appear  to  be  internally 
consistent  when  data  generated  by  a  particular  method  are 
compared . 

Since  the  thrust  of  this  investigation  has  been  the 
comparison  of  techniques  for  predicting  mass-loaded  natu¬ 
ral  frequencies  and  mode  shapes,  and  not  to  argue  the  vali¬ 
dity  of  the  experimental  baseline  data  versus  XASTRAX  data, 
one  would  bo  inclined  to  agree  that  NASTRAX  presents  the 
easiest  technique  for  the  experienced  designer  to  predict 
mass-loaded  modal  data.  The  basis  for  this  is  that  if 
all  the  NASTRAX-predicted  modal  quantities  really  do  exist, 
the  designer  has  at  his  disposal  all  such  data  without  the 
need  for  testing  which  might  load  to  incomplete  results. 

One  area  which  has  not  been  explored  is  the  fact  that 
the  Modal  Analysis  Software  uses  a  least-squares  algorithm 
to  fit  the  Transfer  Functions  for  residues  and  phase  angles 
in  a  given  frequency  range.  In  the  generation  of  the  Modal 
Analysis  Software  mode  shapes,  .my  phase  angles  generated 
due  to  a  frequency  shift  between  Transfer  Functions  were 
ignored.  This,  by  itself,  will  result  in  non-orthoqon.il 
modal  vectors.  Since  MAS TRAN  generates  a  set  of  mutually 
orthogonal  modal  vectors,  it  may  be  possible  to  compare 


the  Modal  Anal vs  is  data  to  the  XASTRAX  data  if  the  effect 


of  phase  angle  (i.e.  -  frequency  shif;  .-/Lth  qrid  point 
within  a  given  mode)  is  included  in  the  display  of  the 
Modal  Analysis  mode  shapes.  This  capability  is  not  pre¬ 
sently  available  on  the  HP5451B  Fourier  Analyser  used  in 
the  experimental  phase  of  this  effort.  Should  this  capa¬ 
bility  be  added  in  the  future,  the  apparent  discrepancy 
between  the  Modal  Analysis  and  MAS TRAM  results  may  bo 
eliminated . 

The  percentage  in  absolute  frequency  shift  from  the 
unloaded  panel  to  the  mass-loaded  data  are  presented  in 
Table  12.  Unlike  Whaley’s  Algorithm  where'  all  but  two 
cases  were  within  +  3°i  of  the  unloaded  panel  results, 

17  out  of  58  (29.3°o)  were  greater  than  this  +  3?£  frequency 
shift.  Even  more  disturbing  is  the  wide  range  of  frequen¬ 
cy  shift  indicated  by  mass-loading  for  some  conf igura tions 
(refer  to  Table  12,  Modes  12,  16,  and  22).  Thus,  it  woul 1 
appear  that  at  least  part  of  the  MASTRAN  data  may  be  ques¬ 
tionable  . 

A  considerable  expenditure  of  the  author's  time  and 
computer  resources  was  made  in  obtaining  these  results. 

For  comparison  purposes,  the  author  kept  records  of  the 
computer  resources  used  to  oblain  the  modal  data  for  each 
of  the  three  finite  element  models  as  well  as  for  designing, 
debugging,  and  running  the  psoudoi nverse  program.  These 
data  are  presented  in  Table  16.  It  should  bo  emphasized 
that  the  author  was  totally  inexperienced  in  the  use  of 
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NAS TRAN  at  the  outset  of  this  effort;  rus,  the  expend- turn 
of  resources  for  the  finite  element  data  using  NAS TRAN 
can  be  expected  to  be  reduced  by  up  to  15%  by  an  engineer 
who  is  more  experienced  in  the  use  of  NASTRAN. 

In  view  of  the  fact  that  the  finite  element  method 
results  did  not  compare  favorably  with  the  experimental  re¬ 
sults,  one  must  answer  the  question  whether  the  time  and 
effort  invested  in  obtaining  those  data  was  commensurate 
with  comparable  investments  from  the  other  two  methods 
vis-a-vis  their  respective  results.  The  answer  lies  in 
the  accuracy  required  by  the  structural  modification,  and 
by  the  vibration  sensitivity  of  the  hardware  to  be  install¬ 
ed.  In  the  case  of  electro-optical  devices  which  are 
known  to  be  inherently  vibration  sensitive,  it  would  ap¬ 
pear  that,  time  permitting,  this  investment  of  resources 
would  be  justified.  However,  one  must  weigh  the  desired/' 
required  results  with  the  method  used  in  each  individual 
case.  Thus,  one  of  the  other  two  methods  (Whaley  Algorithm 
or  Matrix  Recovery)  may,  with  other  less  vibration  sensi¬ 
tive  installations,  be  optimal  when  all  factors  are  con¬ 


sidered. 


Tabic  16.  Computer  Resources  Used 


>tain  Uodt 


Parameter 

Finite  Element 
Models ^ 

Pseudoinver.se 
Matrix 
Recovery ^ 

Central  Pro¬ 
cessor  Unit 
Time  (sec) 

Total 

Average 
Per  Run 

Total 

Average 

Per  Run 

37590 

178.16 

3660 

19.26 

Input/Output 
Time  (sec) 

75155 

356.18 

4640 

24.56 

Computer  Cost 

$4461.46 

$21.14 

$280.92 

$  1.49 

These  data  include  post- processor  time  using  GCSXAST  for 

models  up  to  158  grid  points. 

2  ..... 

The  data  are  for  primarily  nine  grid  point  models  and 

include  post-processor  time  using  DISSPLA. 


The  Use  of  Pseudoinverses  Tn  the  Recovorv  of  the  Discrete 


Matncos,  and  Solution  o 


the  Eigenvalue  Problem 

The  recovery  of  the  mass  and  stiffness  matrices  and 
solution  of  the  eigenvalue  problem  was  considered  by 
Whaley  (Ref  15)  for  a  two  degree-of-f reedom  system  and  th 
case  of  a  square  modal  matrix.  When  one  extends  this  pr 
cedure  to  the  case  of  a  non-square  modal  matrix,  the  pro 
blem  becomes  more  -  implicated  because  the  mass  matrix 
will  in  general  be  'ully  populated  and  thus  the  problem 
becomes  a  generalised  eigenvalue  problem  (Ref  R).  A 
simple  example  for  i  I  iiroc  degroc-of-f reedom  system  is 
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presented  in  Appendix  A.  There  it  is  down  that,  when 
only  two  mode  shapes  and  corresponding  natural  frequencies 
are  considered,  the  resulting  mass-loaded  natural  fre¬ 
quency  and  mode  shape  do  not  closely  approximate  the 
analytical  solution.  Thus,  it  does  not  presently  appear 
that  this  method  will  yield  accurate  results  for  the  case 
of  a  non-square  modal  matrix.  Further  research  into  this 
area  is  needed  as  an  algorithm  (technique)  for  predicting 
mass-loaded  natural  frequencies  and  shapes  for  the  general 
(non-square  modal  matrix)  case  would  be  extremely  useful. 

Summary 

The  time  invested  in  each  method  (including  the  use 
of  Brown's  software.  Ref  2)  was  approximately  equal . 
Although  the  author  does  not  have  such  data  on  Whaley's 
method,  it  is  assumed  since  Glonesk  completed  an  independ¬ 
ent  study  effort  leading  to  a  Master's  degree  that  the 
effort  expended  in  the  use  of  Whaley’s  Algorithm  would  be 
comparable  to  the  other  methods.  Each  method  has  both 
drawbacks  and  good  points,  and  the  use  of  any  one  method 
over  another  will  have  to  be  seasoned  with  a  great  amount 
of  engineering  judgment. 


VII  Recommend' it  i  ons 

In  retrospect,  it  appears  that  the  test  item  selected 
was  much  too  ambitious  for  the  comparison  of  the  three 
modal  analysis  techniques.  One  suggested  test  item  would 
be  a  uniform  thickness  flat  plate.  Many  experiments  have 
been  performed  on  such  an  item  and  analysis  of  the  uni¬ 
form  flat  plate  is  included  in  many  elementary  vibrations 
textbooks  (Ref  5).  Once  the  methods  have  been  tied  to¬ 
gether  using  this  simplified  model,  modifications  to  the 
structure  could  be  added  one  at  a  time  and  the  process 
would  be  repeated.  Thus  one  would  gain  faith  in  the  three 
methods  and,  at  the  same  time,  gain  valuable  experience 
concerning  the  interactions  between  the  various  struc¬ 
tural  components. 

Along  the  same  line  of  the  gradual  increase  in  struc¬ 
tural  complexity  is  the  incremental  increase  in  the  added 
point  masses.  It  is  envisioned  that  this  small  incre¬ 
mental  change  in  added  mass  would  significantly  it  id  in 
tracking  f requoncy  shifts  during  the  method  comparison 
phase  by  allowing  bettor  visual  comparison  of  the  various 
unloaded  and  mass- loaded  mode  shapes. 

The  fact  that  phase  angles  different  from  zero  were 
present  in  the  experimental  data  gives  rise  to  the  ques¬ 
tion  "Were  these  real  modes?"  Speei  f  iru].]y,  this  quest  ion 
is  precipitated  because  a  mode  shape  which  can  bo  visually 
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observed  is  inherently  real-valued  wl'n  no  phase  angle 
between  the  peaks  in  the  Transfer  Functions  at  a  given 
frequency.  If,  in  fact,  a  complex  mode  did  exist,  a  com¬ 
plex  display  would  be  necessary  to  properly  observe  the 
behaviour  of  the  structure.  Further  research  on  the  sub¬ 
ject  of  complex- valued  mode  shapes  is  beyond  the  scope 
of  this  report  but  should  be  pursued  in  the  future. 

One  area  using  the  square/rectangular  modal  matrix 
which  has  not  as  yet  been  investigated  is  the  addition  of 
a  structure  such  as  a  tripod  which  connects  three  or  more 
grid  points.  Hence,  off-diagonal  terms  are  generated  in 
the  mass  and  stiffness  rr.^ilricos  which  further  complicate 
the  issue  because  the  subject  of  finite  element  modelling 
of  the  tripod  structure  (or  any  other  added  structure) 
enters  the  picture.  It  is  the  author’s  opinion  that  this 
would  be  an  interesting  problem  for  one  who  is  interested 
in  furthering  the  research  in  this  area. 
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APPEND 


Mathematical  Considerations  and  Rxarnplo  Problem 
In  the  Calculation  and  Use  of  the  Pscudoinvor  ;c 

Mathemat  ical  Cons  i  derations 

In  the 

pupo 

r  byr  Penrose  (Ref  8)  it  is  shown 

that 

for  any  matrix  A 

there  is  one  and  only  one  matri 

\  A  * 

satisfying 

the  : 

our  conditions: 

(1) 

AA+A  =  A 

(A- la) 

(2) 

AtAA+  =  A+ 

(A- lb) 

(3) 

( AA+  )  *  =  AA+ 

(A-lc) 

(4) 

(A+A) *  =  A+A 

(A- Id) 

The  matrix 

A+,  c 

ailed  the  pseudoinverse  of  A,  is 

the 

matrix  such  that  for  the  inconsistent  sot  of  equations 
A  x  =  b,  the  solution  x  =  A+  B  represents  the  optimal 
solution  to  the  least  squares  problem  A  x  =  b. 

Nov,  consider  the  problem  (posed  by  Noble,  Ref  7) 

A  x  =  b  where  A  =  BC,  and  A,  B,  C  are,  respectively, 
m  x  n,  m  x  k,  and  k  x  n,  and  all  three  matrices  are  of 
rank  k,  then  the  solution  of  A  x  =  b  which  minimines 

a)  the  sum  of  the  squares  of  the  residuals  r  1  r, 
where  r  =  b  -  A  x,  and 

_  t  - 

b)  the  sum  of  the  squares  of  the  unknowns  x  x, 
is  given  by  x  =  A  b,  whore 

A+  =  eT  (C0T)-1  { BTB ) ~ 1  D7.  (A -2) 

When  one  aecompl  Lsie >s  a  docompos i Lion  of  the  matrix  A 
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m 


form  A 


the  form  A  =  IA'  whore  A  is  <ai  m  x  :  matrix,  I.  is 

.  .  .  rP 

an  m  x  m  niatrix,  and  U  is  an  m  x  n  matrix,  B,  B  ,  C, 

rn  rJ'1 

and  C  in  equation  A- 2  can  be  replaced  'ey  I,,  L  ,  U,  and 


mT 


Equation  A- 2  then  becomes  (Ref  10) 


A+  = 

:  UT 

(UUT)  1  (LTL)  1 

1T  (A- 

Example  Problem 

1  . 

Calculation  of  P 

-■pudoi  n verse 

Armed  with 

equation 

A- 3,  lef  us 

consider  an  ex- 

ample  problem. 

Let 

the  matrix  A  be 

represented  by 

1.0 

0.5 

A  = 

2.0 

1.0 

. 

3.0 

0.0 

Following  the  procedure  of  Strang  (Ref  10) ,  the  follow¬ 
ing  procedure  is  used  to  decompose  A  into  L  and  U  factors. 
Using  an  elementary  matrix,  E,  to  accomplish  multipli¬ 
cation/addition  of  rows  in  matrix  A,  and  a  permutation 
matrix,  P,  to  interchange  rows  to  achieve  non-xero  pivot 
elements  in  U  as  needed, 

PEA  =  U,  and  (A- 4) 

A  =  (PE)-1  U  =  LU.  (A- 5) 

The  '  decomposition  of  the  given  matrix  will  now 


be  ealeui  .  . ed . 


- 

- 

- 

r~ 

— 

“ 

1.0 

O.f 

1 

0 

0 

1.0 

0.5 

1.0 

0 . 5 

A 

= 

2.0 

1.0 

- ► 

_  ”) 

1 

0 

2.0 

1.0 

EA= 

0 . 0 

0 . 0 

3.0 

0.0 

-3 

0 

1 

3.0 

0.0 

0.0  - 

1 . 5 

i 

L 

- 

- 

- 

- 

- 

-1 

Not  Lrg 

tha  t 

the  si 

'rond 

row 

of 

EA 

contains 

cl 

n<  n- 

-zero 

pivot  element,  a  row  exchange  with  idle  third  row  of 


EA  is  accomplished  using  a  pormutati  >n  matrix  as 

P  (EA)  = 


1 

0 

0 

0 

0 

1 

0 

1 

0 

L 

1.0  0.5 
0.0  0.0 
0.0  -1.5 


1.0  0.5 
0.0  -1.5 
0.0  0.0 


=  U. 


To  get  L,  augment  the  m  x  m  PE  matrix  vith  an  m  x  m 
identity  matrix  on  the  right  side  and  perform  ele¬ 
mentary  row  operations  until  an  identity  matrix  is  ob¬ 
tained  in  place  of  the  original  position  of  PE .  The 
matrix  now  occupying  the  previous  location  of  the 


original  identity  matrix  is  the  inverse,  (PE) 
[pe!i]  = 


■ 

1 

0 

0 

1 

1 

0 

0 

1 

0 

0  1 

- 

3 

0 

1 

1 

0 

1 

0 

- „ 

0 

0 

i  ; 

- 

2 

1 

0 

t 

l 

0 

0 

1 

0 

1 

0  i 

1 

0 

0 

1 

1 

1 

0 

0 

0 

1 

0 

l 

2 

0 

1 

. 

0 

0 

1 

1 

1 

3 

1 

0 

1 

3 

2 


0 

0 


Thus ,  A  =  LU  becomes 

A  =  LU  = 


r 

r-  i 

r  3 

10  0 

1.0  0 . 5 | 

1.0 

0.5 

2  0  1 

0.0  -1.5= 

2.0 

1.0 

3  10 

o.o  o.o! 

L  -J 

3.0 

0.0 

According  to  Reference  14  if  the  matrix  U  con¬ 
tains  rows  consisting  of  all  zero  elements,  one  can 
delete  those  rows  in  U  and  the  corresponding  columns 
of  L  to  obtain  new  factors  L  and  U  such  that 

A  =  LU  (A- 6) 

Thus , 

A  =  LU  = 


"1 

0" 

[1.0 

0.5" 

1  .0 

0 . 5~ 

2 

3 

0 

1. 

1 0.0 

-1.5. 

2.0 

3.0 

1.0 

0.0 
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The  psoudoinverso  of  A  bocoi^'s 


A+  =  UT  (UUT)  1  (LTL)  1  LT. 


Continuing  to  obtain  A  , 


(A- 7) 


A 


+ 


0.5 

-1.5 


1.0 

0.5 


i.o  o.o 

1.25  -0.75 

14  3" 

12  3 

= 

0.5  -1.5 

-0.75  2.25 

3  1 

0  0  1 

0  0  15 

18  36  -30 


Once  the  matrix  A+  has  been  calculated,  one  must  veri¬ 
fy  that  equations  A-la  and  A- lb  hold.  Since  A  and  A+  contain 
only  real  numbers,  conditions  A-lc  and  A-ld  need  not  bo 
checked.  In  general,  all  four  of  conditions  A-la  through 
A-ld  must  be  checked.  Chocking, 


AA+A  = 


‘1.0 

0.5' 

1 

45 

'  0 

0 

15~ 

1.0 

0 . 5“ 

"1 .0 

0 

.  5 

2.0 

1.0 

18 

36 

-30 

2 . 0 

1.0 

2.0 

1 

.0 

.3.0 

0.0_ 

— 

_  3 . 0 

0.0 

.3.0 

0 

.0. 

^  A. 


S imular lv , 


+  + 

A  AA 


45 


_1_ 

45 


'  0 

0 

15' 

"1.0 

0.5' 

1 

i — ! 

o 

o 

18 

36 

-30 

2.0 

1.0 

1 

18  36  -30 

3.0 

0 . 0. 

*■  i  > 

L  J 

'  0 

0 

15' 

18 

36 

-30 

- 

A  . 

Thus,  the  matrix  A  represents  the  psoudoinverso 


1  he 


original  matrix  A. 

T  + 

If  we  take  B  =  A  and  compute  B  ,  we  find 


+  i  0  18  +  T 

B+  =  0  36  =  (A  )\ 

15  -30 


If  we  let  A  =  U,  B  =  U  ,  equation  3-a  becomes 


U  MU  =  BMA  =  I. 


(A- 8) 


Premultiplication  of  equation  A-8  by  B  and  postmultipli¬ 


cation  by  A  yields 


+  +  +  +  +  + 

B  BMAA  =  B  IA  =  B  A  . 


The  solution  to  equation  A- 9  is 


M  =  B+IA+ 


( A- 9) 


( A- 10) 


provided  the  consistency  condition 


BB+IA+A  =  I 


( A- 11) 


holds.  The  reader  can  verify  equation  A-ll  is  valid. 
Thus,  the  mass  matrix  becomes 


itA*  -  BV  -  ± 


'  0 

18' 

'  0 

0 

15' 

0 

15 

36 

-30 

45 

18 

36 

-30 

~2  648  1296 

43  -540  -1080 


-  540 
-1080 
1125 


'  0.16  0.32  -0.267 

=  0.32  0.64  -0.533 

_-0.267  —0.533  0.556 


( A-l 2 ) 


Before  one  can  be  sure  equation  A-12  represents  the  least 
squares  solution  to  the  mass  matrix,  one  final  check  must 
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be  made; 


B  M  A  -  I .  ('  a-  1  3  ) 

The  reader  can  verify  that  equation  A- 13  holds.  A  similar 
procedure  is  followed  to  calculate  the  unloaded  stiffness  and 
damping  matrices  using  equations  2  and  3,  respectively. 

Example  Problem  2.  Calculation  of  Mass-Loaded  Mass  and 
Stiffness  Matrices 

Following  the  procedure  of  Example  Problem  A-l,  consider 
the  three  degree-of-f reedom  system  of  Figure  A-l.  Let  Mi  = 

M2  =  M3  =  1  and  Hi  =  K2  =  K3  =  100.  Neglecting  damping  and 
the  forcing  functions,  the  equations  of  motion  become: 


~1 

0 

o' 

/  •  • 

*1 

'  200 

-100 

0  ' 

N  0 

0 

1 

0 

*2 

+ 

-100 

200 

-100 

U  =  p 

0 

0 

1 

(x3 

0 

-100 

100 

(X3I  p 

These  equations  are  the  three  degree-of-f reedom  eigenvalue 
problem  for  this  system.  Upon  solving  equation  A-ld  for  its 
eigenvalues  (natural  frequencies)  and  eigenvectors  (mode 
shapes),  the  results  in  column  two  of  Table  A-l  are  obtained. 


Table  A-l.  Analytical  Data  From  Throe  Dogreo-of-Freedorn 
Spring-Mass- Damper- System 


Parameter 

Mi=M2=M3=1 

k1=k2'k3~100 

M1=M3=l,  M2-2 

Ki=K2=K3= 100 

Unnormalized 
Modal  Matrix 

m 

r,j  2=12.47 
f  jl)  3  —  18 .019 

iip3.813 

0)0  =  11  •  8'lrj 

0)3=15.6-18 
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f  1  ( t ) 


X 

x3 

5SS 

■ 

S m 

M2 

\1 

"3 

'-l 

i 

m 

'-2 

bid 

'-3 

QZ2 

Equations  of  Motion: 


Mass  Matrix  Damping  Matrix 


"Mi  0  0 

Hi 

(CX+C 2  “C2  0 

(*l) 

0  m2  0 

x7 

+ 

-C2  (C0+C3  -c3 

x2 

0  0  M3 

*3 

0  -c3  c3 

|x3l 

Stiffness  Matrix  Forcing 


(Ki+K2 ) 

-k2 

0 

1*1 

fl(i 

~K2 

(Kd+Kd ) 

-k3 

1  x2 

= 

f2 

0 

-k3 

K3  J 

( x3 

f3(t) 

Function 


Figure  A-l.  Three  Degree  of  Freedom  System 


Next,  let  a  mass  of  magnitude  one  be  added  to  mass  M2  sucl 
that  M2  is  now  equal  to  two.  This  will  simulate  a  mass- 
loaded  configuration.  Data  for  this  configuration  arc  listed 
in  column  three  of  Table  A-l.  Note  that  the  mass  addition 
lowered  the  natural  frequencies  and  modified  the  mode  shapes. 
Now,  consider  the  first  two  natural  frequencies  and  cor¬ 
responding  mode  shapes  (Table  A-l,  column  2)  to  be  the  meas¬ 
ured  modal  data  in  the  frequency  range  of  interest.  Using 
the  procedure  of  Example  Problem  A-l  to  recover  the  mass;  and 
stiffness  matrices,  the  "generalized  inverse"  equations  of 
motion  become: 


0 


7 


r 


"  .485  -  .  105  -.131" 

xl| 

j  J.  —  •  :  .  _  1  2  1.  .  /  1  9 

1 1 

-.105  .353  -.237 

*2  + 

-34.212  47.009  -22.523 

x. 

-.131  -.237  .249 

•  •  1 

x3| 

_  11 .719  -22.523  12  .  S71_ 

x3l 

i°l 

0  .  (  a-  1  5 ) 

l°l 

After  adding  a  unit  of  mass  to  the  M22  term,  and  resolving 
the  eigenvalue  problem,  the  following  results  were  obtain¬ 
ed: 

W1  =  0.02 
(jJ2  =  4.45 
=  12.46 

Comparison  of  these  results  with  column  two  of  Table  A-l 
reveals  that,  with  the  exception  of  o' ^ ,  the  values  are 
nearly  identical  with  the  unloaded  results,  and  that  the 
magnitude  of  the  frequency  shift  was  nowhere  near  that 
which  was  encountered  in  column  three  of  Table  A-l. 
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APPENDIX  B 


APPENDIX  n 


Detailed  Modal  Analysis  and  Test  Procedures 

Modal  Assurance  Criteria  (MAC)  Funct ion  Testing 

The  test  panel  was  suspended  from  a  ceiling-mount¬ 
ed  conduit  in  the  Structural  Vibrations  Branch  (Bldg  2 -1C) 
of  the  Air  Force  Wright  Aeronautical  Laboratories  (FBG/ 
AFWAL)  by  bungy  cords  which  were  approximately  four  feet 
in  length.  These  bungy  cords  were  attached  to  the  panel 
through  holes  in  the  skin  located  near  each  of  the  four 
edges.  The  test  grid  used  by  Glonosk  (Ref  3)  was  marked 
on  the  upper  surface,  and  a  reference  accelerometer  was 
mounted  on  the  skin  (see  "x"  on  Fig  2)  using  double-backed 
tape.  The  reference  accelerometer  and  moveable  accel¬ 
erometer  responses  were  individually  amplified  and  filt¬ 
ered  through  a  pair  of  matched  filters  before  their  signals 
were  input  to  the  HP-5451B  Fourier  Analyser  (Ref  ■!). 

The  reference  accelerometer  was  used  to  provide  the  (rig¬ 
ger  signal  for  data  sampling. 

With  the  panel  configured  as  described,  t  biO  move¬ 
able  accelerometer  was  mounted  at  grid  point  number  one. 

The  panel  was  then  excited  (tapped)  13  times  at  randomly 
selected  locations.  Following  each  tap,  and  prior  to 
computations  using  these  data  in  the  Fourier  Analyser, 
each  data  sample  was  viewed  separately  on  a  ref resh-t  r.«re 
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oscilloscope*  to  insure  the  tapping  h  d  :  ot  overloaded 
the  internal  electronics,  or  that  multiple  excitations 
had  not  occurred  during  the  data  sample.  Assuming  an 
acceptable  data  sample  had  been  obtained,  these  data 
were  stored  in  the  Fourier  Analyser.  Each  time  a  good 
data  sample  was  processed,  the  response  data  for  that 
sample  was  averaged  with  the  previous  samplo(s)  such  that 
after  15  iterations  the  stable  average  for  each  channel 
was  computed.  Using  these  average  response  data  with 
the  reference  accelerometer  as  the  input  and  the  move¬ 
able  accelerometer  as  the  system  forced  response,  the 
Fourier  Analyser  calculated  both  the  MAC  and  Transfer 
Functions  for  that  grid  point  and  displayed  the  Transfer 
Function  on  the  oscilloscope.  At  this  point  the  user 
had  the  option  to  either  accept  the  data,  or  reject  the 
data  and  reaccomplish  this  process  for  that  grid  point. 

If  accepted,  hard  copies,  including  plots  and  diciitiul 
iniormation,  were  made  of  Transfer  Function  and  MAC 
Function  data.  The  Transfer  Function  data  were  stored 
for  later  use. 

The  next  stop  was  to  successively  mount  the  move- 
able  accelerometer  at  each  of  the  remaining  grid  points 
and  repeat  this  procedure.  Once  the  entire  grid  had 
been  surveyed,  the  MAC  Function  data  were  examined  'o 
determine  candidate  f requeue ies  for  whi ch  a  cu r vo  fit 
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of  the  Transfer  Function  data  would  be  :  x:e.  Prior  to 
curve-fitting  the  Transfer  Function  data,  each  discrete 
Transfer  Function  was  examined  at  the  given  frequency 
range  to  determine  the  data  with  the  "cleanest"  spike. 

This  spike  was  fit  first  using  a  least-squares  algorithm 
to  determine  the  natural  frequency  and  damping  ratio  to 
be  used  in  each  of  remaining  Transfer  Functions  to  deter¬ 
mine  the  mode  shape  at  that  frequency.  This  process  was 
repeated  for  each  candidate  frequency  range1  identified  by 
by  the  MAC  Function  data  in  the  range  of  0-500  Hz. 

Modal  Analysis  Test  Procedure 

Several  grid  sets  were  used  in  this  testing:  the 
25  grid  point  scheme  of  Glenesk  (Ref  3,  Figure  2),  a  25  grid 
point  scheme  (Figure  5)  to  coincide  with  the  addition  of 
structural  components,  and  a  65  grid  point  model  with  fine 
meshes  superimposed  on  each  of  the  four  internal  bays  (Figure 
6)  to  investigate  the  modes  of  vibration  of  each  of  these 
bays.  Since  the  tost  procedure  is  identical  for  each  qrid 
set,  it  will  only  be  discussed  for  the  25  grid  point  s'Mtvho. 

This  testing  used  a  moveable  accelerometer  to  survey 
the  qrid.  However,  unlike  the  MAC  fund  ion  testing,  a 
second  stationary  accelerometer  w  is  not  used.  Instead, 
one  grid  point  was  selected  (number  one  for  this  grid) 
at  which  the  structure  would  bo  excited  by  an  impact 
hammer  configured  with  a  load  coll  to  measure  t  ho  amplitude 
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of  the  impact.  The  output  from  t !,.  >  It,  ;  .  roll  was  used  a.-, 
the  trigger  source  with  the  moveable  accelerometer  meas¬ 
uring  the  forced  system  response.  Again,  unlike  the  ‘-'.AC 
function  testing,  the  structure  was  excited  only  five  time 
at  each  grid  point.  The  number  five  was  chosen  for  two 
reasons:  first,  it  was  noticed  during  the  MAC  function 

testing  that  more  than  five  data  samples  had  no  signifi¬ 
cant  impact  on  increasing  the  accuracy  of  the  average  I’SD; 
and  second,  due  to  the  large  number  of  grid  points  used  in 
the  larger  models,  including  more  than  five  iterations  per 
grid  point  would  cause  excessively  long  testing  while  pro¬ 
viding  little,  if  any,  improvement  in  the  quality  of  the 
data . 

The  process  began  by  mounting  the  accelerometer  at 
grid  point  number  one  and  exciting  the  structure  using 
the  impact  hammer  five  times  at  the  selected  excitation 
point.  As  in  the  MAC  function  testing,  the  data  wore 
examined  following  each  discrete  sample  to  determine 
if  an  overload  of  the  internal  electronics  or  a  multiple 
excitation  had  occurred.  Assuming  acceptable  data  wore 
obtained,  the  discrete  force  input  was  multiplied  a 
force  window  (Figure  B-l)  to  artificially  force  it  to 
die  out  rapidly  as  would  be  trio  case  for  an  impulse  ex¬ 
citation.  similar]},  the  structural  forced  response 
was  multiplied  by  an  exponent,  i  ally  deriving  window  (!  iquie 

B-2)  to  simulate  the  effect  of  damping  on  the  structural 
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The  average  of  cosh  s  t.mp.L<  •  ..  !  h  proceeding 
samples  was  computed  after  each  of  the  five  discrete  dut.  : 
samples  was  accepted.  As  before ,  u  Transfer  Function  was 
computed  for  each  point.  Upon  completion  of  the  fifth 
iteration,  the  average  Transfer  Function  and  a  Coherence 
Function  wore  computed  and  individually  displayed.  A 
decision  whether  to  accept  the  data  was  based  on  a  Coher¬ 
ence  Function  with  many  of  its  values  at  or  near  one  in 
the  0-500  Hz  range.  This  process  was  repeated  until  the 
entire  grid  had  been  surveyed. 

Once  the  survey  was  complete,  a  grid  point  was  select¬ 
ed  which  was  believed  to  contain  all  the  modes  of  the 
structure.  Several  methods  for  examining  the  Transfer 
Function  data  at  that  grid  point  were  available  in  the 
software.  Among  *  hese  methods  were  a  Kennedy-Plancu  Circle, 
fit,  a  least-squares  algorithm,  and  an  option  which  merely 
computed  the  magnitude  of  the  Transfer  Function  at  -electee! 
f req  enrios  where  the  real  and  imaginary  parts  of  the  Trans¬ 
fer  Func*. ion  wore  90  degrees  out  of  phase  with  each  other. 
Initially  the  author  selected  the  magnitude  option  as  it 
ignores  damping.  However,  the  .least-squares  a  Igor  it  hr.  v.  as 
even  t  u  a  1 1 y  u  s e  d  to  a  n  a 1 ys  e  the  d  u i a . 


Ma  t  hem  a  t  i  eg  1  -  F  xpo  r  i  mont  a  1  i’roce  dure 

This  procedure  user,  the  experimentally  do f  e rrr: i  ned 


modal  d.'.t  a  for  the  unloaded  panel  gathered  i  ri  the  Hod  il 
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Analysis  Section  presented  earlier  in  this  Appendix.  The 
mathematical  portion  of  this  procedure  begins  with  the 
general  matrix-vector  differential  equation  used  in  modal 
analysis , 

[M]  x  +  [C]  x  +  [Kj  x  =  f(t),  (B-l ) 

where  [M]  is  the  mass  matrix,  [C]  is  the  damping  matrix, 
and  [K]  is  the  stiffness  matrix,  all  of  which  are  square 
and  symmetric  matrices.  A  generalized  coordinates  solu¬ 
tion  to  this  equation  results  in  "n"  uncoupled  equations. 
Assuming  the  solution  to  equation  1  can  be  written  in  the 
form 

x(t)  =  [U]  q(t)  ,  ( B-2 ) 

where  [U]  is  a  square  matrix  consisting  of  mode  shape 
vectors  as  its  columns  (i.e.-the  modal  matrix),  then  it 
follows  that  q(t)  is  the  generalized  coordinate.  Following 
the  derivation  of  Meirovich  (Ref  5)  after  substitution  of 
equation  B-2  into  equation  B-l  and  premultiplication  of 
equation  B-l  by  [U]T,  the  following  definitions  are  made: 


[U]T[MJ  [U]  =  [I] 

Tx  °" 

[U]T[K]  [u]  =1 

[u]T[c]  [u]  =|o  2?n(,)n 


(B-3a) 


( B-3b) 


(B-3c) 
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With  these  definitions  equation  B-l  becomes 


qr  ***  2  f  r^r^r  +  -x'r^5r  =  fr*r=:lj2,3>...  ( B— 4 ) 

m 

where  fr(t)  =  [U]  f(t)  is  the  generalized  force. 


Equation  B-4  can  be  solved  for  at  most  "n"  natural 
frequencies.  However,  there  are  often  less  than  "n" 
natural  frequencies  in  the  frequency  range  of  interest 
such  that  the  modal  matrix  contains  fewer  columns  than  it 
contains  rows.  It  is  proposed  to  premultiply  each  of  the 
equations  B-3  by  ^[U]T^  and  postmultiply  each  of  these 
equations  by  [u]-1  in  such  a  way  as  to  isolate  the  mass, 
damping,  and  stiffness  matrices  on  the  left  side  of  equa¬ 
tions  B-3  which  involve  these  respective  matrices.  In 
general,  [U]T  and  [U]  will  be  non-square  (rectangular) 
matrices  necessitating  the  use  of  the  pseudoinverse  devel¬ 
oped  by  Penrose  (Ref  8)  to  find  their  inverses.  Using 
a  superscript  "+"  to  indicate  a  pseudoinverse  and  after 
the  previously  mentioned  pre-  and  post-multiplications, 
equations  B-3  become 


[M]  [U]  [U]+ 
[K]  [U]  [U]+ 

[C]  [U]  [U]+ 


^[U]T^  [I]  [U]+  ( B-5a ) 

([U:|T)  („'<!)  [u]+  (B"5b) 
(tU]T)  (q2^11  )  tu:!+  <B‘5c)' 
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Equations  B-5  can  be  further  reduced  using  standard  linear 
algebra  matrix  manipulations  (Noble,  Ref  7)  to  yield  the 
mass,  stiffness,  and  damping  matrices  respectively  as 

[M]  =  ([U]T)+  [I]  [U]+  (B-6a) 

tc]  -  (Cu]T)+(o"2f 

Reference  8  (Penrose)  contains  the  basic  theory  of  the 
pseudoinverse . 

With  these  definitions  for  the  mass,  damping,  and  stiff¬ 
ness  matrices,  and  the  experimentally  determined  modal  data 
for  the  unloaded  panel  gathered  in  the  Modal  Analysis  Section 
presented  earlier  in  this  report,  one  can  calculate  the 
mass,  damping,  and  stiffness  matrices.  Appendix  A  contains 
a  simple  example  of  the  technique. 

Since  the  thrust  of  this  procedure  is  to  accurately 
predict  the  mass-loaded  natural  frequencies  and  mode  shapes 
of  a  complex  structure  when  only  the  unloaded  modal  data 
are  known,  the  procedure  must  be  continued  to  recover  mass- 
loaded  modal  data  n>  Un).  For  example,  consider  the 

mass  mafrix.  When  the  designer  has  made  his  decision  re¬ 
garding  the  placement  of  the  additional  mass  (hardware)  in 
the  structure  in  question,  he  need  only  add  an  appropriate 


n  n 


") 


[U]  + 


(B-6c)  . 


°)  [u]+ 


( B-6b) 
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mass  in  the  mass  matrix  at  the  proper  grid  location  to 
obtain  the  modified  mass  matrix.  A  similar  procedure  is 
followed  to  obtain  the  modified  stiffness  and  damping  ma¬ 
trices.  Assuming  one  is  only  interested  in  mode  shapes 
and  natural  frequencies,  the  modified  mass  and  stiffness 
matrices  are  substituted  into  equation  4.  Equation  4  is 
then  solved  for  the  mass-loaded  natural  frequencies  and  mode 
shapes . 


Finite  Element  Modelling  Procedure 

Finite  Element  Models.  The  three  finite  element  models 
used  in  this  analysis  were  described  in  Section  II.  Because 
the  development  of  a  finite  element  model  is  described  in 
Reference  11,  the  main  question  became  how  to  add  a  quanti¬ 
ty  of  mass  to  the  mass  matrix  which  would  be  used  to  solve 
for  the  modified  natural  frequencies  and  mode  shapes. 

NASTRAN  provides  this  capability  in  the  form  of  a  C0NM2 
card  which  adds  a  finite-valued  point  mass,  including  its 
inertia  properties  and  neutral  axes  offset,  to  a  particular 
location  in  the  mass  matrix.  Another  card,  the  ASET1  card, 
was  used  to  discard  all  but  the  out-of-plane  translation 
from  the  analysis  set  in  the  solution  for  mode  shapes  and 
natural  frequencies.  Appendix  C  contains  listings  of  a  sam¬ 
ple  NASTRAN  deck  used  in  the  modal  analysis  of  each  con¬ 
figuration. 

Analytical  Procedures.  NASTRAN  (Ref  11)  uses  several 


rigid  formats  to  provide  flexibility  to  the  user  in  the 
analysis  at  hand.  Rigid  Format  3  was  chosen  for  the  modal 
analysis  of  the  test  panel.  This  format  neglects  damping 
and  solves  equation  B-l  with  the  forcing  function  and  damp¬ 
ing  matrices  set  equal  to  zero.  The  output  from  this  For¬ 
mat  is  a  table  of  the  "n"  natural  frequencies  and  tables  of 
the  "n"  eigenvectors  {mode  shapes).  These  data  were  then 
post-processed  using  GCSNAST  to  be  displayed  on  a  standard 
computer  terminal. 
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FB'.TlOO,  l  01 00  *C Ml **0000.  T 7O0YT4, AY* |  NSCN.^o*; 


4TTACH,|MSL ,10-LlBRAR Y ,SN-ASD. 
iTTACM.rmSRLA  ,  IP-L  IBRARV,SN-4SO. 

LIBRARY,!  -SL  .0  IS  SPL  A. 

FTN.PL-IOOOO. 

LCO. 

ATTACM,PLT103B,ID“A780?6T,SN"ASD. 

RLT1018. 

RE  Ml  NO,T 4PF9*I, 

ROUTE ,TAPFB9,TID-AF , F I D-F B l , S T- CS A , DC  - P R . 

"  ENO  OF  RECORD 

PROGRAM  matrix* input, output , rl FILE-0*TAPE  *- INPUT, TAPE  6-OUTPUT  > 

c* 

CP 

CP 

CP  THIS  IS  a  TEST  VERSION  OF  MATRIX  USING  THE  NINE  GRID  POINT  MODEL  ANO  THE 
C*  RESULTING  "CDf  SHAPES  ANO  NATURAL  FREQUENCIES  TO  VERIFY  THE  VALIOITY  OF  THF 
C*  PROGRAM  TO  CALCULATE  NASS-LOAOFO  NATURAL  F  RfGUF  NC  If  S  ANO  *»OOE  SHAPES.  NINE 
C*  SEPARATE  masS-LOADED  CASES  ARE  CONSIOfRFO  IN  ACCDROANCF  WITH  THF  CASES 
C*  TESTFO  BY  GLfNFSK.  SEVERAL  CASfS  FOR  EACH  CONFIGURATION  ARE  CONSIDERED: 

C*  THAT  IS,  SEVERAL  "STARTING  POINTS"  APE  CONSIDERED  FOR  THE  RANGE  OF  N A  T  UR  — 
C*  At  FREQUENCIES  TO  SIMULATE  A  REAL-WORLD  CASE  WHERE  ONE  IS  ONLY  INTERESTED 
CP  IN  A  SPECIFIC  FREQUENCY  RANGE. 

C* 

C* 

C* 


c* 

C • TH I S  PROGRAM  SOLVES  THE  CENFRAL  T R I x-VE C TOR  DIFFERENTIAL  EQUATION 
CPFOR  THE  PASS,  DAMPING,  ANO  STIFFNESS  "ATR ICES  USING  GENERALIZED 
CPCOOROINATES.  THE  FOLLOWING  OEE1NITIONS  ARE  * AOE : 

CP 

C*  UT  •  AN  •  U  •  AI 

C*  UT  P  AC  P  U  -  ZETA 

C*  UT  •  AK  P  U  •  ONEGA 

C* 


c« 

WHERE  U 

THE  POO  AL  NATRIX 

c» 

UT 

■ 

THE  TRANSPOSE  OF  THE  NODAL  «ATR!X 

c» 

AN 

■ 

THE  NASS  MATRIX 

c* 

AI 

m 

THE  IDENTITY  MATRIX 

c* 

AC 

m 

THE  OAPPTNC  MATRIX 

c« 

AK 

■ 

THE  STIFFNESS  NATRIX 

c» 

ZETA 

m 

A  DIAGONAL  * A  T  R I X  WITH  THF  PRODUCT  OF  THE 

c* 

NATURAL  FREOUfNCY  ANQ  THE  DAPPING  RATIO  FOR 

c» 

EACH  NOOE  OROEBED  IN  INCREASING  PODE  NUMBER 

c* 

ALONG  THE  N  A  I N  DIAGONAL  AND  ZFROES  ELSFWHERE 

c* 

ONEGA 

m 

A  DIAGONAL  MATRIX  WITH  ThF  SQUARE  DF  THt 

c» 

NATURAL  FREQUENCY  FOB  EACH  PODf  ALONG  THE 

c« 

*AIN  DIAGONAL  ORDERED  IN  INCREASING  NODE 

c* 

NUMBER  AND  ZEROES  ELSEWHERE. 

CP 

CP 

CPSOLVING  THFSF  EQUATIONS  USING  STANOARD  ■  AT  R  !  X  *•  A  N  I  PUL  A  T  I  ON  S  YIELDS: 


CP 

CP 

AN 

•  UGITN 

p  Al 

P  IJGIN 

CP 

AC 

•  UGITN 

♦  ZETA  P  UGlN 

CP 

AK 

•  UGITN 

♦  P* 

'EGA  ♦  tIC,  1 N 

CP 

CP 

WMfRF  UG|N 

-  THF  GENERALIZED 

INVERSE  OF  THF  NORMALIZED 

CP 

POOAL  MATRIX 

CP 

UGITN 

•  THE  TRANSPOSE  OF 

UGlN 

c* 

CPNOTE 

TO 

the  user: 

IE 

J  IS  THF  NUMBER  OF 

GRID  POINTS  AND  L  IS 

THE 

CP 

nunrer  of  «onrs  in  the 

FREQUENCY  RANGE  OF 

INTEREST 

CP 

THE 

"ATRICFS  OFF  INFO  *B0VF  SHdlRD  BF  D  1 BF  NS  I  ONF  D 

CP 

AS 

FOLLOWS: 

CP 

AN,AC*AK,AI  ,ZETA,  ANO  ONEGA:  DIMENSION 

J  X  J 

CP 

UGlN 

:  dimension 

L  X  1 

CP 

UGITN 

:  OlNENSION 

1  x  L 

CP 

U 

t  DIMENSION 

1  X  L 

CP 

ZCTA5,OPFGAS,AL 

:  DIMENSION 

L 

CP 

CP 


Figure  C-l.  Program  Matrix  Listing, 
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r. 


C*N#0  •  NUMERICAL  ROW  DIMENSION 
C*NR  •  NUWMFR  OF  RQWS 
C*HC  •  NU-NFR  CF  COLUMNS 
CP|NlTIALt/F  THf  MATRICES: 

. . . 

D I  “F N$ TON  UUf9,9t,0!FI9t ,01Ff |9) ,01F If  91 ,01E?(91 

01  Mf  NS  ION  XA(3t,YAU!,7Af3f,X(9f,Y(9t,  70(91, XX(3l,Yr  (31, ??n> 
DIMENSION  XXX(9|,YYY!9I,77M9| 

Ot-FNSION  A  Al P  I  9  | 

0l«fNS  19*  UM9,9>  ,U6(9, 9|,U7|9,9>  ,U8f9  ,9) 

DIMENSION  WM*0A|9| , M ALP  19 I  , U A  I 9  ,  9  I  • F R R OR ! 9 I 

DIMENSION  AA  n  9,9  I  ,F(9,9|  ,019*91  ,H(9,9  I  ,0NECA(9,9I  ,U(9,9J,uU9,9| 
DIMENSION  ACt9,9l,Al(9,9|,7ETA(9,9t,FF(9,9l,Gr.t9,q|,HM|9,9I.MA(9,9 
U,FFE  (9,9)  ,GGGt  9,9),  C(9, 91  ,0(9,91,8(9,91 , 7 E T A S I  9 > ,0 ME C AS 1 9  I , ALl 9  I 

dimension  erporu9) 

DIMENSION  NAMB0A|9|,UC1TN|9,91 
OMENSION  SU*M9|  ,SUMH9»,SU"?I9»,SUM3<9t 
DI-ENSION  SUMAI9I 

REAL  AM«fq,9l  ,U3P(0,9|,S1<M,UK3(  16  I « U3 ( 9  ,* Y , RE T A? I  9 t , WK * ( 18? I 
REAL  A(9,9f,UGINf9,9),T0l ,S|9) ,WM !62),AMf 9,9) ,AK|9,9I ,9tTAt9) ,U2I 
19,9),V;K7U*7I  ,AAN|9,9|  .AAKI9.9)  ,RET  All  9)  f  MK  1  1  1671  ,U1  19,9)  ,W"EGA(9  ) 
CO -RLE*  ALFA  191 ,719,9)  ,ALF  A1I9  1, 7U9,9|,0MfnA0<9|,BL  AM0OAI9I 
COMPLEX  ALEA?19) ,7319,9), ALPI9I 

DIMENSION  ONEGlf  8, 81 ,OmFCAII 8),A1I8,8),CCI9,8I,U3TP|9,8I,AMA|9,9I, 
10019,8) 

NR -9 
NC»9 
ND*M 

CALL  C0MPRS 
NCASES-3 

DO  N  RK-1,NCASES 

. . . 

c* 

C*  READ  STRUCTURAL  GRID  POINTS 

CP 

DO  20?  I  —  1 ,  NR 

REAO«,Y(I l,Y(I 1,70111 

fFf  FOFISLINPUD.NE.Ol  STO® 

20?  PRINT*, *Mt,YM), 70(1) 

. . . . . 

CP 

CP  INITIAL I7E  THE  MATRICES 

CP 

C*P*P*P*PPP****'**#*P'PPPPPPPPPPPP«PPPPPPPPP**tPPPPPP«PPPPPPPPPPPPPPPPPP 

DO  111  l-i,NR 
00  111  1-1, NC 
All, J) -0.0 
Ul  I, J)»0.0 
U2M.  Jt-0.0 
UA  II,  JI-0.0 
UT(I,I)-0.0 
UG1NI J,I 1-0.0 
UG I  TNI  I, JI-0.0 
2ETASI J)-0.0 
ALIJI-l.O 

111  OMEGAS! J)-0.0 
00  2  l -l , NC 
00  2  J-l.NC 
AMI,  JI-C.O 
OMEGA t I «  J1-0.0 
.?  2ETAII  ,  Jl-0.0 
00  333  1-ltNR 
00  333  !-l,NR 
AMf I, JI-0.0 
AMI,  JI-0.0 
333  ACM,  JI-0.0 
00  70  1-1, NO 
00  70  J— I , NO 
AIM.  JI-0.0 

to  o*  rcu i,  Ji-o.o 

PR  TNT*  ,  — THE  INITIAL!  7F0  MATRICES  ARE:** 

PRtNTP,  "  THE  MOOAL  MATRIX:- 
DO  M  1-|, NR 

81  PRtNTP, fU! I • J), J-1,NC) 


Figure  C-l.  Program  Matrix  Listing  (Continued) 


POINT**"  IMF  "GOAL  "  A  T  0  I  x  TRANSPOSE:" 

no  67  J ■  1  * NC 

67  POINT*. IU1I J, II. 1-1. NO* 

POINT*,"  THf  CfNFOAL 1 1  f  0  INVERSE:" 

cn  n  i-i. nc 

6J  POINT*. IUCINIl.il. 1*1. NOT 

POINT*,"  THF  GENEOALI7EO  INTI  OFF  TO  ANSPOS  E  :  " 
00  64  1-1. NO 

64  PRINT*, IUGITNII. 11. 1-l.Nfl 

POINT*,"  THE  FETA  NATOlx:" 

00  6S  1-1. NC 

65  PRINT»,IZt TAII , II , 1*1. NCI 

PRINT*,"  THF  0"FGA  NATO IX:" 

00  bb  1-1. NC 

66  POINT*. I0NFGAI1, II. I-l. NCI 

PRINT*,"  THf  «OU"NX*  «ATO|Cf S:" 

POINT*,"  The  AL  "AT01X:" 

on  in  i-i. nc 

111  PRINT*,A1II1 

POINT*,"  THE  **  6S  "IATOIX:" 

00  117  1-1. NO 

117  PR  I  NT*. I *■« I . 11 . 1-1 .NOI 

POINT*,"  the  OA"P|NG  HATO|X:" 

00  111  I-l. NO 

113  POINT* ,f ACI I . II . J*1 <NR 1 

PRINT*."  THE  STIFFNESS  "ATOIX:" 

00  11*  I-l. NO 

114  PRINT*,! AK I  I  , 1 1 •  1- 1  .NO t 

POINT*,"  TUI  Il'TNTITY  "ATOIX:" 

00  11S  t-l.NC 

115  PR  [NT*.  I  AM  I  ,  II.  1- I  ,NCt 
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P  R I N  T • , SU - f  Jl 
SU-I  JI-SQRT  (SU*<  Jl  I 
00  ?\ *  1-1, NR 
Z  18  AC  I,  11-41  I,  jl/MIRf  Jl 
00  7?l  I-l.NR 

ZZl  PRINT#, (A||, J|, J-i, NC I 
00  ??0  I-l.NR 
00  ZZ 0  J-l,NC 
01 1,  ll-RI I • Jl 
22 0  UTI  J,l  l-A( I, Jl 

c* 

C  ♦  Ofl  f  TF  THf  U5T  CPIU-N  FRO**  THf  MODAL  MATRIX  U  TO  f  OK  -  THf  RfCTANGUlAR 
C*  -ATRIX  in. 

c* 

c««# ••♦•••*•••♦♦•♦♦♦••♦•♦♦♦•♦•♦•♦••♦♦•••••♦*♦••♦*♦ ••♦*•♦♦♦♦••♦•♦♦•♦•♦♦•♦♦♦•• 

00  81  I - 1 tNR 

on  fli  j-i.no 

81  UT (f*JI-A(|«J) 

PRINT*, -THIS  IS  THf  U  -ATRJXt- 

00  NU  I  ■  1  ,  NR 

911  PR  TNT*  ,(UI  I  « Jl  * J«1 «NR  I 

PRINT*,-TH!S  IS  THf  A  HATRIX:- 
00  91  I-l.NR 

91  PRINT*, (All, JI,J-1,NR| 

PRINT*, -THIS  IS  THf  UT  MATRIX:" 

00  9?  I-l.NR 

92  PRINT*, (UTII, Jl, J-1,NR  I 
PRINT*, -THIS  IS  THf  U1  MATRIX:" 

00  8 Z  I-l.NR 

8  2  PRlNTP.fini!  ,  Jl  ,  J-l  ,N0I 

00  7 1  1*1, NO 
71  AL  P ( I  1-0,0 

on  81  i-i, no 

0-FG1 t 1,1 l-ONf CMC I »♦•* 

83  Aid,  11*1.0 

PRINT*, -THIS  IS  THf  0*FG1  -ATRIX:" 

00  84  I -  1 , NO 

84  PRINT*. (P8F  GUI,  Jl,  J-l, NO! 

PRINT*, -THIS  IS  THf  A  1  MATRIX:- 
00  8S  I  —  1  ,  NO 

85  PRINT*, 1  AH  I  ,  ||  ,  J-l, NOT 

(;••«#•**•**«*#••#**♦•***♦*•*•****•*•**•*•*♦****+•*********•*♦*«***«•••••••** 

CP 

CP 

CPCONPUTf  THF  GFNFRAU/FO  INVFRSf  OF  THF  MODAL  HAT  R I X  ANO  ITS  TRANSPOSE 

CP 

CP 

C*P****#***P*P*«* **••**•• 

NR01-9 
NR  0-9 
NC  *9 
NR  -9 
TOl-O.O 

?0  CALI  LGINF(A,NRO,NR,Nr,Tni.,UCIN,NRPl,S,)<K,IfRI 
PRINT*, -THIS  IS  THE  MATRIX  UGlNt- 
00  77  1-lfNC 

77  PRINT*, l UG  tNI|,J),J-l,NRI 
00  173S  l-l.NP 
00  1714  J-l.NC 
173S  UC1TNI I , JI-UGfNI  J, I  I 

PRINT*, -THIS  IS  THf  MATRIX  UGI TN  :" 

DO  414  I-l.NR 

414  PR  f  NT  #  ,  I IIC  ITNfl.JI.J-t.NCI 

PRINT*, "THE  VAlUfS  (IF  7E  T  A  ANp  O-fCA  A  R  f ,  RESPECTIVELY:" 

00  IS  f-1, NC 

08FGASC  T I -0-f CAS  I  It / l op. 

IS  PRINT*,/fTAS(  | l,P-4  CAST  II 

. . . . . ***** . . 

CP 

CP  CONPUTF  THf  GFNFRAVWTO  INVERSE  OF  THf  -ATRIX  til.  THIS  IS  THf  MATRIX  UTP. 

CP 

C**#********P*«*P**«*********** *«*•***•**•*♦*•«•••• . •♦••*♦•♦#••••♦•♦•• 

NR  01 -9 
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NR  3-<» 

NC3-* 

70L-0.0 

NRD4-8 

CALL  LCINF |U3 ,NR0  3,NQ  3 ,NC3,T0L  »U3P,NR04.S  3,WX  3, If  R  3 ) 
PRINT#, "THIS  is  1HF  "A*9!X  U 3P ? " 

00  30  l  •  1  »  NO 

30  PR 1NT#, IU3PC ! , it  •  J-l ,NC» 

c# 

C#  TRANSPOSE  tl3  TO  FPR9  U1T.  ALSO  TRANSPOSE  U3P  TO  FOR*  U  3  T  P. 

C# 

00  31  I -l, NR 
00  31  l-l,N0 
3  t  U3TP | I , M-U3P(  J, f \ 

PRINT#, "THIS  is  IMF  -ATRJY  u  3  T  P  :  - 
DO  3?  I  •  1  » NR 

3?  PRINT#, IU3TPI I , Jl , J-l, HOI 


C#FOR"  THE  7E  T  A  ,  OPEGA,  ANO  Af  MATRICES? 

C# 

C# 

(#############♦♦##*##'###########•######•••#' 

PR  I  NT  #, "THE  ACTUAL  "ATRICFS  ARE:* 

00  5  l-|,NC 

*eTA(I,f l-?#2f TASI I I#0»tcaS I  I > 

PEGAI I, l  I -0#E  GAS  f I »##2 
Af (I,  1 1-1.0 

5  cont'nue 

PRINT#,-  THE  7E  T  A  p  ATr  f  X ;  — 

DC  ?5  l-l.NC 

75  PRINT#,  (  7f  TAM,  M  I  •  I  I  •  l  ,NCt 

PRINT#,-  THF  0#E  C  A  PATRIY:- 

00  ?6  1-1, NC 

26  PRINT#,  (OPf  CAM  ,  1 1  >,!!•?, Nt> 

PRINT#,-  THE  IDENTITY  -ATR  Its* 

00  IT  I  -  I  ,  NC 

27  PRINT#, I A1 1 1, Jl , J-l, NCI 
00  4i  l-t.NR 

00  41  1 T-l ,NC 
CM  ,111-0.0 
DO  41  R-1,NC 

CM,Hl-CM,m«UCITNM,Kt#AIfK,lI) 

41  CON  r  fNUf 

00  42  I-l,NR 
00  42  If  - l , NR 
AN  M , l I  1-0.0 
00  42  K-1,NC 

AM!  ,1  II  -  A  PM  •  1 1  MCI  I  |-N1#UG1N<  K,  1 1 » 
AA«(I , ! I  )  -  A* I  I  ,1  I  1 
4?  CONTINUE 

00  43  1-1 , NR 
00  43  !!-l,NC 

oiifin-r.o 

00  43  *-l,NC 

Ofl*m-riI,IIMUCITNM,RI#7FTAIK,m 

43  CON  T INUF 

00  44  I • 1 , NR 
00  44  I  I • l , NR 
ACM  ,1  M  -0.0 
00  44  K-l.NC 

AClItlll-ACII •! 1 l#Oll ,K»#UG1N|K, | l» 

44  CONTINUE 

00  45  I  - 1 , NR 
00  45  tl-l.NC 

eu*m-o.o 

00  45  K-l.NC 

6M,TM-EM,Ml#UCITNM,KI#0PeCAIK,  IM 

45  CONTINUE 

00  46  1-IfNR 
DO  46  t I • I , NR 
AMItl  I  *  -0.0 
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UC ITNI l  .K l*UC I N( K  ,11 


00  *6  K-l.NC 

rkii.iii-pkii.iii-eu.ki-uginik.iii 
46  CONTINUE 

00  1  I-l.NB 
00  l  l-l.NC 

IFt8SSI8Ml.JM.LE.10.E-8l  I'll,  11*0.0 
If I  485  I 4K I  I « 11 I  •  Lf • IO.f-8)  INI  1 ,  11-0.0 
1  If  IARS  UC  I  I  •  1 1  I  ,L£  .  10.E-8  I  »C(I.  11*0.0 
PRINT*.”  THE  D»«PINC 

00  48  1-1. NS 

48  PRINT*,IACI1,1I»1”1*NRI 

PRINT*, ”  THE  STIFFNESS  NlTRI*:” 

00  40  1-1, NR 

49  PMNT*’.”CHSCK  ^(INSISTENCY  CONOITION  FOR  IHF  10ENTITY  H4TRI*: 
00  000  I  - 1 .NR 

00  000  1-1. NR 
f I  I, 11-0.0 
00  000  K-l.NC 
000  f  tl,  I! -MI.1I 
00  00?  1-1. NR 
DO  00?  l-l.NC 
HI  I ,  I  1-0.0 
00  OO?  K-l.NR 

00?  HI  I ,  1 1 -HI  I .  1 1  ♦  f I  I  .K I  *UIK ,  1  I 
00  001  l-l.NC 
on  oo3  i-i.NC 

4  41 11 •  II -0.0 
on  001  K-l.NR 

003  *»l I  I . II ■**! 1  I » I*  ♦  UT I  I ,K I *H( K , 1 1 

00  0003  l-l.NC 
00  0003  l-l.NC 

0003  IF  <  4  85  1 4  4  I I  I  ,11  I.LE.  10.E~8  I  4A11I.JI-0.0 
00  004  l-l.NC 

004  PRINT*. I  4*11  I .  II •  1*1 .NCI 

PRINT*. "CHECK  CONSISTENCY  CONOITION  FOR  THE  ?ET4  N4TRU- 
00  90S  l-l.NC 
00  OOS  l-l.NC 
Ff II,  11-0.0 
DO  OOS  K-l.NR 

005  ff I  I ,  ll-FF I  I  .11  ♦  UTII  ,KI*uC1TNIK,1I 
DO  006  l-l.NC 
00  006  l-l.NC 
CGI  I  .  11 *0.0 
or  006  K-l.NC 

006  CCI  I  .  11 -CCI  I  .  i)  *  f  F  1  I  ,K  i*;etmk  ,  II 
OP  00?  l-l.NC 
OC  00  7  l-l.NC 
HHI I. 11-0.0 
00  00 7  K-l.NR 

007  HHI  I ,  I  I -HHI  I  , 1  I  *  UCINI  I  ,K I *U« K •  II 
00  008  l-l.NC 
00  008  l-l.NC 
Hi  I  I ,  II -0.0 
00  008  K-l.NC 

008  H4  I  l  ,  II -H4I  I  ,  II  ♦  CC 1  I ,KI «HMI K . II 
00  0004  l-l.NC 
00  0004  l-l.NC 


0004  If  U8SIH1II.1I1.1E.I0.E-8I  H4II.J1-0.0 
DO  000  l-l.NC 

000  PR INT«, IHAI I , II ,  I- l.NC I 

PRINT* , "CHECK  CONSISTENCY  CONOITION  fne  THE  08EC4  «ATRIX:« 
DO  1003  l-l.NC 
00  1003  l-l.NC 
FFf I  I ,  11-0.0 
00  1001  K-l.NC 

1003  fff f I,  ll-f f f 1  I .11  ♦  f F  I  I.KI*nNEC*«K.II 
00  1001  l-l.NC 
00  1001  l-l.NC 
CCCI I .  II -0.0 
00  1001  K-l.NC 

1001  CCCI I ,  I  l-CCCI  I . II  ♦  M f I  I ,K l*HH|K .1 1 
00  0005  l-l.NC 
00  0005  l-l.NC 

0005  If lARStGGG 1 1 . II  I .C F . 10.E-P I  CCCII, 11-0.0 
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- tra  "m'rM.hC - - 

100*  PRINT*.  fCCCK.  Jl. J-l, NCI 
PR  I  NT ♦, "T  HE  NASS  NATRix:* 

00  100*  1*1. NR 

100*  PRINT*. (  ANU  .  Jlf  J-l, NCI 

C*********************** «*#•#*********•*♦*»*»***»»*•«♦••»•*•»*••♦»'•* •"•*•**• 

c* 

c* 

c* 

C*  SINCE  THE  NASS  NAS  ADCEO  TO  GRID  POINT  NUMBER  FIVE  ON  THE  ACTUAL  PANEL.  THE 
C*  NASS  N  A  T  R I  X  WILL  BE  PERTURBED  IN  THE  S.S  POSITION  R>  ADDING  A  QUANTITY  Of 
C*  NASS  EQUAL  TO  THE  NASS  LOADING  ON  THE  ACTUAL  PANEL.  THIS  KILL  SJnuLATE  THE 
C*  "REAL  WORLD"  CASE.  NOTE  THAT  THE  NASS  WAS  PLACED  ON  ONLY  ONE  GRID  POINT.  If 
C*  THE  NASS  WAS  A  STRUCTURE  SUCH  AS  A  TRIPOD  WHICH  COwNrCUO  TO  THREE  GRID 
C*  POINTS.  THE  ADDED  STRUCTURE  WOULD  HA VF  TO  BE  "ODFLLFO  IN  THAT  IT  WOULD  AOD 
C*  BOTH  NASS  ANO  STlffNESS  TO  THE  PANEL.  THIS  ADDED  NASS  AND  STIFFNESS  nOULO 
C*  GENERATE  ADDITIONS  TO  Off  DIAGONAL  TERNS  OF  LIKE  SU^PfR  TO  THE  TwO  CONNECT- 
C*  EO  CRIO  POINTS.  AS  PREVIOUSLY  NOTEO,  THIS  PROGRAN  SOLVES  THE  GENERALl/EO 
C*  EIGENV At UE  PROBLEN.  THUS.  THERE  IS  NO  NEED  TO  CREATE  A  S 1 NUL AR  -OOEL  FOR 
C*  ANY  ADDITIONAL  OANPING  INCURRED  BY  THE  ADDITIONAL  STRUCTURE. 

C* 

C* 

PR INT* ."ENTER  THE  LOOP  TO  PERTURB  THE  NASS  NATPIXt" 

IfCKW.FQ.ll  ANfS,Sf-AAMS,SI«G.ll** 

IF  IKK. EQ. 2  I  AN(  S.  S  I-  A A* (S.S  MO.?  *20 
lfIKK.fO.3l  AN(SvS»-4AW(S,S|«0.*lSe) 

PRINT*."  TMf  NASS  naTRIx:" 

DO  SO  I-l.NR 

SO  PR  INT*.  I  AAM  I  ,  n  ,  J-l  ,NC» 

PRINT*. "THE  NOOIflfD  NASS  "  A  T  R  I  X  :  * 

DO  *7  I-l.NR 

*7  PR INT*. ( A"l 1 •  Jl  .  J-l . N« » 

PRINT*. "THE  STlffNESS  "ATR1X:" 

DO  SB  1-1. NR 

SB  PR INT*. C AKI I , J I, J-l.NC I 

C************** ****** «*•«*•**•*********•*****•**•**•*•»****••**•**•*• ********* 
C* 

C*  SOLVE  THE  GENERALIZED  EIGENVALUE  PROBLEM  fQR  THE  NASS-LOAOfO  NATURAL  FR£- 
C*  OUENCIES  and  "OOE  shapes  with  The  Square  nodal  NATRIX  u. 

c* 

C************************* ••**•******♦*♦*♦•*♦**<*♦♦*♦*♦♦♦****♦***•**•*••••• 

1AN-Q 
IAK-9 
NANAK-9 
12-9 
I  JOB-2 

CALL  E  IG2FIAK, IAK.AN.I AN.NANAK .1  JOB. ALFA. BETA. 7*1/ »WK2*l£R2! 

N-9 

P«INT*.*THFSE  ARE  THE  VALUES  OF  RLANBOA:" 

00  123*  I-l.N 

RLANBDA I  II -CSQR T I AL F A  I  I  I /BETA! II) 

PRINT*. RLAweOAC I  I 
123*  CONTINUE 

PRINT*. "THESE  ARE  THE  VALUES  OF  THE  2  BAiRIVS* 

00  1009  I-l.N 
1009  PRINT*.C7(I. J). J-l.N) 

PRINT*. "THE  PREOICTEO  "ODE  SHAPES  ANO  NATURAL  FREQUENCIES  FOR  ANIS 
l.SI-".A*(S,SI ."ARE :• 

00  21  t-l.NC 

21  NANBOAI  I  l-CSQRTIRL  ANBOMI  »**2I 
PR INT* .IWANBOAf I )» l-l .NCI 
LA-9 

CALL  VSRTAIWANBOA.LAI 
PRINT* .IN A"BOA(l I, l-l, NC) 

00  22  t-l.NC 
00  22  J-I.NC 

IFIABS IWAWPOAI I l-RFALI RLA"RDAf  J)  II  ,LE.  1O.E-0 )  GO  TO  23 
IFlABSfWAPBOAl ll-REALC RL  AN80A I J| I I.GT.10.E-B)  GO  TO  22 
23  00  2*  K-l.NR 

U2CK, I l-C.O 

2*  U?«X.  !  t-U?<K,1l*RE ALC/CK, Jll 

22  CONTINUE 

00  IS?  J-l.NC 
SU«*I JI-0.0 
00  1  S3  I-l.NR 

IS 3  SUN*!  JI-SU-M  JT  *U?I  1  .  Jl*»2 
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SU«Ht  Jl-SQRTfSUHM  Jit 
00  15?  1-1, HR 

152  U2tl .JI-U2C I, JI/SU«4< i» 

PRINT*, -COMPUTE  SU«  OF  SOS  OF  DIFFERENCES  BE  T  WE  E  N  LOADED  ACTUAL  AN 
10  LOADED  PPFOICrCD  MODAL  VfCTORS:- 

. . 

c* 

C*  READ  mass-LOAOEO  MODAL  VECTORS 

C* 

00  250  |-1 ,NC 
DO  ?50  I  •  1 ,  NR 
UUll, JI-0.0 
RE  A0l5,251  I  UUI I , il 

251  FQRMATtlOX, €15.51 
250  PRINT#, UUI !« Jl 

00  25?  I  •  1  ,  NR 
00  25?  1-1, NC 

252  UUII, JI-UUI I, J»*UUI I, JI/100. 

00  263  J ■ 1 ,NC 

SU"l tit-0.0 
00  265  1-1, NR 

265  SUNK  Jt-SUPll  JI*UUt!  ,  Jt##2 
SUM! ( Jl-SORTf  SUM l { j|| 

DO  263  1-1, NR 

263  UUC r ,  Jl-UUtf , il/SUNl I  Jl 

c«, 

c* 

C*  COMPUTE  THE  SUM  OF  THE  SQUARES  OF  ThE  DIFFERENCES  BETWEEN  ACTUAL  AND  PRE- 
C*  OICTED  LOAOFD  MODAL  VECTORS 

C* 

. . . 

00  253  1-1, NC 
OIF  t  Jl-O.o 
00  253  |-1, NR 

253  OIF  I  J I -01  Ft  JtMUUt  t,  J)-U?tf  ,  it  I*  *2 
00  259  I-l.NC 

OIFltlt-O.O 

260  OIFlt 1 l-DIFI II 
259  CONTINUE 

00  261  1-1, NR 

261  OIFim-SQRTf  OIFlt  I  »» 

PRINT*, -THESE  A«E  THE  VALUES  OF  THE  MODE  iHA PE  VECTORS:- 
00  ?55  1-1, NR 

255  PRINT#, |U2C I, Jl • J-l,NCt 
00  86  I • 1 , NR 
00  86  J- 1 ,N0 
CCtt, JI-0.0 
00  86  K-1,ND 

86  CCt  I,  Jt-CCt  I,  JMU3TP  I!  ,K1*Alf K,  J» 

DO  87  I - 1 , NR 

DO  87  1-1, NC 
AMAU,  Jl-0.0 
ANMf I , JI-0.0 
00  87  K-1,N0 

AMMtl  ,  Jl-AMWt I , J|*CCI I  ,KI*U3PtK,  Jl 

87  AHAtl, ll-AMNt 1 ,it 
00  88  I  - 1 , NR 

00  88  J- l , NO 
OOt  I, JI-0.0 
on  88  K-1,ND 

88  OOt I , JI-OOI I • Jl *U3TPf I ,RI AQMEClf K , J | 

on  89  I  - 1 , NR 

00  89  1-1, NC 
AAM  !•  JI-0.0 
00  89  K-l.NO 

89  AAKtl, J I -A  AM  II  fJt«COII  .RI*U3P|R,JI 
00  72  I-1,NR 

00  72  1-1 , NC 

Iff  ANSI A--t I , Jl f.Lf . 10.F-8I  amm||vji-o.O 

1F1  A8SI  A-A(  I  ,  J  I  1  .  L  F . 1 0 . E -5 •  A-M  I,  Jl-0.0 
72  IFJABStAAM!  ,J»  J.lF  .iO.f-ftt  AAKII,  JI-0.0 

C***** . . . . 

C* 

C#  THIS  IS  THE  IFAST-SOUARFS  SOLUTION  TO  THE  MASS  MATRIX  BEFORE  THE  NASS  HAS 
C*  BEEN  AOOEO  Tn  THF  5,5  Tfa*. 

. . . 

Figure  Ol.  Program  Matrix  Listing  (Continued), 


PH  INT  *,"TH|  S  IS  THE  «AT8|X  *M: 
00  7 3  l-l.NH 

PR ! NT •»  t  *•'"1  1 «  ft  ,  J  - 1  ,NC  I 
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c** 

PRINT*. "THIS  IS  THE  NATR  | X  AAK:" 

00  78  I  *1  *  NR 

7*  PRINT*,! AAKtI. Jl • J-l.NCI 

. . . 

C* 

C*  SOLVE  THE  GE  Nf  R  AL  I  7  E  0  EIGENVALUE  PROfllEN  F  UR  T  HE  "ASS-lOADEO  NATURAL  FRE- 
C*  OufNClCS  ANO  HOOF  SHAPES  USING  THE  RECTANGULAR  NODAL  MATRIX  U3. 

C* 

()•♦***•#*♦***♦**♦**•***•••***•***•****•*****••*•***•••••**•**••*••••*•****•*• 
If  IKK.  EQ.  il  ANN |  8,81-  ANAI  8,8) *0. 1 18* 

IF  IKK.EQ.71  ANNI8,8I-Ana|S,8I*0.78? 

IF  IKK. FQ. 3)  ANN|5,St-ANAIS,SM0.*lS8 

17-9 

I  JOH-7 

IA-9 

IB-9 

N-9 

PRINT*, "THE  NOOIfICD  NASS  "ATRIX  AN«:» 

DO  6?  I-l.NR 

67  PR|NT*,I APN|| , Jt , J- l ,NC I 

CALL  EIG7FIAAK,|A,ANN,IH,N,Ijnp,AiFA7,liETA7,71*I7'Wt'8,IER8l 

C****************** ************** *••***••*•••**•«*****•****•*****•••*•• ****** 

c* 

C*  THESE  ARE  THE  IE  AS  T-SQUARE  S  SOLUTIONS  FOR  THE  NATURAL  FREQUENCIES  WITH  THE 
C*  THE  AOOFO  NASS  IN  THt  8.5  POSITION  OF  THE  NASS  NATRIX. 

C* 

c********* ♦**•••♦•♦•••♦• ••••••*•••*♦*•••••♦*•♦*•♦•♦••♦•♦♦♦♦•♦*•••••♦•♦♦•••••♦ 

PRINT*, "THESE  ARE  THE  VALUES  OF  ThE  73  NATRIX:" 

00  7*6  I-l.NR 

7*6  PRINT*, 1731  I. Jl, J-l.NCI 

PRINT*. "THESE  ARE  THE  VALUES  OF  ALP:" 

00  13  I-l.NR 

ALP  I  I  I -C SORT f ALf A7  ! 11 /8FTA?< I M 
WALPU  T-CSQRT  «ALP  (I»**7I 
31  PRINT*. ALP!  I  I  , WAl P I  I  ) 

CAIL  VSRTAINALP.LAI 

PR  INT*, "THESE  ARE  THf  VALUES  Of  kJHP:m 

PRINT*, IWALPl I T.I-l.NCI 

00  765  1-1, NC 

IFI  I . E  0. 1  I  AALPI9I-WALPI I  I 

IF(I.GT.1I  GO  TO  766 

GO  TO  765 

766  K-I-l 
AALPIKI-WALPI  I  • 

765  PR 1NT*,AALP(KT 
DO  198  I-I  ,  R 
00  198  J-l.NC 

IF  f  APSlAAlPd  l-RFAL  (ALP<  JIII.LE  .1O.E-0I  GO  TO  197 
tF IANS! AALPI I l-REAL I ALPI Jl » I.GT.  10. E-8  I  GO  TO  198 

197  00  196  K-l.NR 

UA|K,1  i-o.o 

196  U  8  I  K  ,  I  I  •  IJ8  f  K  ,  I  I  »  R  E  A  L  f  7  1  f  K  ,  J)  I 

198  PR INT* ,  U8 I  I .  I  I 
00  180  I  ■  1  . NC 
SUN1I JI-C.O 

00  181  l-i.NR 

181  SUN1I  JI-SU-H  JI»U8I  I  ,  J»**7 
SUN)|  JI-SORTf  SUNK  J|  | 

00  180  I-l.NR 

180  U8l  I  ,  IT  -U8I  I  •  J  I /SUNK  f  • 

PR INT* ,"CC"PUTf  THE  SUN  RF  SQS  OF  DIFFERENCES  HE  T  W  M  N  LOADED  ACfUA 
IL  ANO  PREOICTEO  LOAOFD  NODAL  VECTORS  USING  THE  GENERALIZED  INVERSE 
7:" 

00  788  l-l.NC 
OIFFI JI-0.0 
00  788  I-l.NR 

788  OlFMJT-niFft  JI-IUUI  I  ,  JI-U8I  1,  Jl  |**7 
DO  786  I-l ,NC 
787  0  IF  7 ( I 1-riF E I  I • 

786  CONTINUE 

00  767  l-l ,NR 

767  DIF7I f l-SGRTIOI F7C 111 

.  Proqram  Matrix  Listinq  ( ‘  ’ont  inuod ) 
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PRINT*, "THE  EIGENVECTORS  ARF:" 
on  ?8  1-1. NR 
PRINT*, (U2(l  .  J).  J-l  .NCI 

PRINT*,*TME  GENERALIZED  INVERSE  EIGENVECTORS  ARE:" 

00  79  I  •  l  ,  NR 

PRINT*. t  U*  t 1 • 1 1 . J-l.NC  I 

00  8*  I-l ,NC 

A  At.  P  I  II-  A  AL  P  (  1 1  *100. 

NAHBOAl ! | -HANBO AC  I  1*100. 
on  52  L8-1.NC 
REAOC5.53I  WM8DAIL01 
FORNATI27X.F8.3l 
PRINT*. W«BDAILfl» 

00  55  !»1 ,NC 

ERROR! !l-( IWNBOA  f | l-WANOOAt I  1 1 /M"bDA(!  I  I *100. 

00  231  1-l.NC 

ERROR  1  (n-<(wNftDA(It-AALPI!n/NM8DAUl  1*100. 


C*  COMPUTE  THE  OOT  PROOuCT  OF  THE  ACTUAL  MASS-LOADED  VfRSUS  P»FDICTEO  NASS- 
C*  LOAOEO  MOOAL  VECTORS  FOR  THE  SQUARE  MOOAl  MATRIX  !U7I  AnO  FOR  THE  RECTAN- 
C*  ULAR  MOOAL  MATRIX  tUBI 

C* 


00  1*5  I-l, NR 
00  1*5  J-l.NC 
05 1 ( , J! -0.0 
U6f I • Jt -0.0 

U5U  .  JI-UM  I.  JI*U2<  J.  I  I 
1*5  U6I1  .JI-UM  I.  J»*UU(  J.  I  I 
00  1*6  I-l, NR 

on  i*6  j-i,nc 

U7I I . JI-0.0 

unit# JI-0.0 

DO  1*6  K-l.NC 

U7(  I  ,  JI-U7U  .  jl  «UUU  ,K»*UMK,  J  I 
1*6  U8I I. Jl-uei 1 • J»*UU< I ,K  »*U6CK , J| 

DO  1*7  I-l  ,  NR 

1*7  PR  INTI6, 1*4.81  107(1,11.  J-l, NCI 
UU  I**'*  I-l, NR 

1**9  PR  fNf  (6, 1**8  I  <UM  1  ,  II  ,  J-l  ,NCI 
1**8  FORMAT! 2X,0f 2X,F 7. 3 || 

PRINT*, "THE  EIGENVALUES  ARE:" 

PRINT*, "THE  COMPARISON  OF  F I GENVALUES : * 

PRINT*,"  UNtOAOEO  PREDICTED  ACTUAL 

IRE  PREDICTED  P  E R C F N T  SOUARF" 

PRINT*,"  EXPERIMENTAL  "ASS-LOAOFO  maSS-LOAOED 
IT  MASS-LOADED  ERROR  ROOT" 

PRINT*,"  DATA  USING  EXPERIMENTAI 

l  USING  PREOTCTFO  OF" 

PRINT*,"  SQUARE  MOOAL  DATA 

1*5  PSEUOOINVERSE  VERSUS  OIFF’S" 

PRINT*,"  MATRIX 

10  C9R0WS.8C0LSI  ACTUAL  S Q  *0 " 


PERCENT 


EXPERIMENTAL  PREDICTED 


PRINT*,"  SQUARE  "GOAL  DATA  VERSUS  OIFF 

1*5  PSEUDOINVERSE  VERSUS  0!FF*S" 

PRINT*,*  *  A  TR  IX  ACTUAL  SQ* 

10  C9R0WS,8C0LSI  ACTUAL  SQ'O" 

DO  29  l-lfNC 

OMEGAS! 1 1 -OMEGAS  II 1*100. 

PRINT! 6,51  I  O-EGASn I. WAMPOA ( I  I , MMfiQA  f  II.FRPORlII.DIFlII  I  . AALP (  I  I  , 
1ERR0R1 f I » ,0IE2 ( I  I 

FORMAT!*X,F7.3,7X,F7.3,7X,f 7.3,6X,E6.2 ,*X,f 5.2»5X,F7. 3. 7X.F6.2  ,*X, 
IF 5 •  2  I 

TFILLL.FO.il  GO  TO  269 
IFfLLL.F0.3l  GO  TO  2*9 


C*  IF  LLL  F0UAL5  ONE,  THIS  STEP  WILL  RE  I CNQfi F D  ANO  Thf  PROGRAM  WILL  CO  TO 
C*  STATEMENT  26 9,  HOWEVER,  IF  LLL  EQUALS  ?,  THE  wasS-lDADFO  *004^  VECTORS 
C*  I  ACTUAL!  WILL  BE  READ  IN  AND  THEIR  MQDE  SHAPES  PLOTTED. 

C* 


2S8  00  268  J-l , NC 

00  268  I  - 1  , NR 
REAni5,267l  U2fl.il 

267  F0RMAT!13X,E:0.2I 

268  U2(t, JI-U2II, JI/100. 
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DO  1*0  J-1«NC 
5U*?t JI-0.0 
00  1*1  1-1.  N* 

1*1  SUK?  1 J  I  -SU*21 Jl  *U2  1  I. J I  ••? 

SU*2! JI-SQRTl$U*2l J> I 
00  1*0  1-1. NR 

1*0  U2 I! « JI-U2I I « J1/SUN21 II 
269  CONTINUE 

DO  211  LL-l.NR 
C *LL  BCNPLIH 

CALL  TlTtim-wOOE  .  C  ONF  l  CUR  ft  T 1  ON  S* , 100 . 3 . 0 , 8 . 0 1 

CALL  AXFS301M*ETA3",lOO,"7%*,,lOO,"RAOl  AL3*,l00,l5.  ,13.  ,  13,  1 

CALL  VUAASI-20.* 70.. SO. > 

CALL  G*AF 301-20. .10. *50. , -1 .  , 10.. 50. *-0.5 .10. .SO. 1 
LK-Ll 

CALL  NFMRfNIll 
00  21?  KL-1,2 
1F1KL.F0.1I  CO  TO  ?0* 
lFlKL.CT.il  CO  TO  201 
20*  DO  206  1-1.3 
K-3MI-II  ♦! 

206  CALL  CU0V3PU  |K|  ,Y1K1  W01KI  ,3.1  I 
00  207  J-1,3 

K- J 

KXll I-K1K I 
YT1 1  I-Y1KI 
f  1 1 1 1- IQfK I 
X*f2l-Xf  KOI 
YYf  2>-Yf  KOI 
22  f?  •  -  201  KOI 
XX 1 31-Y|K«6l 
Y V | 3 | -Y 1 K ♦ 6 1 
2213I-701K46I 

PRINT*#XX|l!fYYll|,77m,XXl?|,YY12l,Z?121,XXf  3I.YY13I  ,22  1  31 

207  CALL  CUR V3 01 XX, YY, 22,3,11 
IFfKL.FO.il  GO  TO  203 

201  00  203  L-l.NR 

XXX1LI-X1L) 

YYY1LI-YIL1 

IFfLLL.FO.ll  CO  TO  199 
IFILLL.F0.2I  CO  TO  199 
1FILLL.F0.3I  2/7UW0<Ll*UA(L«LLt*5. 

CO  TO  203 

199  72211  WOft»*U2fL,LL»*3. 

203  PRINT*, XXXfLI, YYYfLl  tZZMLI 
205  DO  209  1-1,3 
K-3M  1-1 1*1 

209  CALL  CURV301XXXfKl,YYYfKl,2?ZfKl,3,ll 
00  210  J-1,3 

K»  ) 

XMlI-XXXfKJ 
YAf 1 1-YYYfK | 

7ai  1 1-  n  mk  » 

XAI2I-XXX1K43I 
YA  f  2  I  -  YV  Y  f  K  ♦  3  | 

Jh\ ?l-272fKO» 

XA f 31-XXXf K*6l 
YAf  3I-YVY1K46I 
ZA|3l-7mX*6l 

PR  1  NT*  ,X  i  1 1  )  ,  YA  f  1 1 , 7  A  f  1 1  ,  XA  f  2  I  t  Y  A  1 2  I  v  2  A  f  2 1 «  XA  f  3  I  ,  Y  A  f  3  I  ,  Z  A I  3  I 

210  CALL  CORY  301XA.YA.7A ,3,11 
fFIKL.CT.il  CO  TO  216 

208  CALL  NOHIOF 
CALL  NfWPFNflf 
UfKL.FO.ll  GO  TO  217 

216  CALL  NFWPFNU1 

217  CONTINUF 
CALL  OASH 

212  CONTINUF 

211  CALL  FN0PLI1I 
6  CONTINUF 

CALL  OONFPL 

213  STOP 
FNO 

1  FNO  OF  *FCf?®0 


Figure  C-l.  Frogicitn  Matrix  Listing  (Continued). 
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-10.  24.3 74  12.  5 
-5.5  24.174  12.5 

10.  24.374  10.4 
-10.  6.  12.4 
-4.4  6.  12.4 
lO.  fe.  10.4 
-10.  0.  12.4 
-4.4  0.  12.4 

' 

I  10.  0.  10.4 

l  44 

l 

1*,6.617  7 

.1589694 

5 

1 

2  00 

2 

181.781  7 

.2807813 

4 

7 

1  103 

0 

706.  !*»«. 

. 7570244 

5 

1 

4  101 

0 

706.66  3 

.711 3849 

5 

4 

2  110 

0 

2  36.0  30 

.6190009 

5 

2 

4  no 

0 

716.717 

.  7678094 

5 

5 

4  tio 

1 

270.003 

.7539196 

5 

4 

2  22  4 

1 

451.660 

.7763402 

5 

2 

3  2  42 

1 

686.667 

.7087252 

5 

3 

1  1 

168. 

.OOOf *00 

.OOOf  *00 

30.5 

•  0 

.0 

1  2 

7.60 

.OOOf  *00 

.OOOF *00 

72.6 

•  0 

•  0 

1  1 

-24.6 

.ooof *oo 

.OOOE  *00 

94.5 

.0 

.0 

1  4 

81.3 

•000E *00 

•OOOF *00 

63.  1 

•  0 

•  0 

1  4 

-44.2 

.OOOf *00 

•  OOOf  *00 

69.5 

.0 

.0 

1  6 

-31.2 

.OOOf *00 

.OOOf *00 

22.  7 

.0 

•  0 

1  7 

10.4 

.OOOf  *00 

.ooof *00 

108.9 

•  0 

•  0 

1  6 

4.30 

• OOOf ♦ 00 

.QOOf  *00 

54.1 

.0 

.0 

l  o 

-11.3 

• OOOf ♦ 00 

.ooof  *00 

90.  7 

.0 

.0 

2  l 

132. 

•  OOOf  *00 

•000E*00 

121.4 

•  0 

•  0 

2  2 

2.41 

.ooof  *oo 

.OOOF *00 

.  3 

•  0 

,  0 

2  1 

-27.6 

.ooof  *oo 

.OOOF  *00 

165.7 

.0 

•  0 

2  4 

127. 

.OOOf *00 

.OOOf *00 

105.  1 

•  0 

•  0 

2  4 

-45.0 

.OOOf *00 

.ooof  *00 

122.5 

•  0 

•  0 

2  6 

-23.4 

.OOOf  *00 

.OOOF  *00 

28.  1 

.0 

.0 

2  7 

4.41 

.ooof *00 

.ooof *00 

172.1 

.0 

•  0 

2  8 

8.30 

.OOOf *00 

.OOOt  *00 

4.  3 

.0 

.0 

2  0 

-10.5 

.OOOF *00 

•  OOOf  *00 

139.2 

.0 

.0 

l  1 

80.8 

.ooof  *00 

.OOOf  *00 

177.1 

.0 

.0 

1  2 

-14.6 

.OOOf  *00 

• ooof ♦CO 

122.1 

•  0 

•  0 

l  3 

151. 

.ooof  *00 

.OOOE  *00 

171.6 

.0 

.0 

1  4 

700. 

.OOOf *00 

•  OOOE  *00 

141.6 

.0 

.0 

l  4 

202. 

.ooof  *oo 

.OOOf  *00 

140.4 

.0 

.0 

l  6 

-780. 

•ooof *oo 

.OOOf *00 

135.8 

.0 

•  0 

l  7 

-142. 

.ooof *00 

.ooof  *00 

117.1 

•  0 

•  0 

1  8 

-24.1 

.ooof *00 

.000? *00 

170.1 

•  0 

•  0 

1  o 

143. 

.ooof *oo 

,000€  *00 

120.6 

•  0 

.  0 

4  1 

48.0 

.ooof *oo 

,000f*00 

45.2 

.0 

.0 

4  2 

-11. 0 

.OOOf *00 

.OOOE  *00 

111.6 

•  0 

.0 

1*8. 

.ooof *oo 

•OOOE ♦ 00 

108. S 

.0 

•  0 

4  4 

737. 

.OOOf *00 

.0006*00 

114.', 

•  0 

•  0 

4  4 

228. 

.ooor*oo 

■  o  oi)  f  ♦  0  0 

131.6 

.0 

.0 

4  6 

-275. 

.ooof  *oo 

.ooof *00 

126.  7 

.0 

.0 

4  7 

-121. 

.ooof *oo 

•000E*00 

122.6 

.0 

•  0 

4  8 

-25.4 

.OOOf *00 

. ooof ♦oo 

119.2 

•  0 

.0 

4  0 

140. 

.OOOf *00 

■  OOOE  *00 

11  3.4 

.0 

.0 

2  1 

264. 

•ooof *oo 

.OOOE *00 

64.  5 

•  0 

.0 

2  2 

-14.0 

.OOOf *00 

•  OOOE  *00 

150.2 

.0 

•  0 

2  1 

142. 

.ooof *oo 

•  OOOF  *00 

141.9 

.0 

.0 

2  4 

77.2 

.ooof *oo 

•  OOOE  *00 

98.5 

.0 

.0 

2  4 

-131. 

•  OOOf  *00 

•  OOOE  *00 

145.9 

.0 

.0 

2  6 

-402. 

.ooof *oo 

•  OOOE  *00 

140.0 

.0 

•  0 

2  7 

-06.4 

.ooof *oo 

•  OOOE  *00 

148.  7 

.0 

.0 

2  6 

-3n.  8 

.OOOf *00 

.OOOF  *00 

149.  7 

.0 

•  0 

2  0 

1  12. 

.OOOf *00 

.OOOE  *00 

142.6 

.0 

.0 

4  1 

146. 

.OOOF  *00 

.OOOE *00 

56.8 

.0 

.0 

4  2 

8.51 

.ooof*oo 

.OOOF  *00 

29.0 

•  0 

•  0 

4  3 

87.7 

.oocr  *oo 

.OOOF  *00 

161.5 

.0 

.0 

4  4 

11.6 

.OOOF  *00 

.ooci *00 

96.  7 

•  0 

.0 

4  4 

-101  . 

•  OOOF  *00 

.OOOF  *00 

160.6 

.0 

.0 

4  6 

-417. 

.0001  *00 

.OOOF ♦no 

*  17.  6 

.0 

.0 

4  7 

61.6 

.000*  *00 

•  OOOF  *00 

1  1.4 

.0 

.0 

4  8 

-29.6 

.  100f  *00 

•  OOOF  *00 

153.8 

.0 

•  0 

4  0 

81.  1 

.OOOF *00 

.ooof  *oo 

160.  \ 

.  0 

•  0 

4  1 

421. 

.OOOF *00 

.ooof  *00 

22.  1 

.0 

•  0 

4  2 

-11.0 

• OCOF *00 

.OOOF  *00 

123.9 

.0 

.0 

%  1 

124. 

•  OOOF  *00 

•  OOOF  *00 

l  74.4 

.0 

.0 

4  4 

80.8 

.OOOF  *00 

•  OOOF  *00 

21.1 

.0 

.0 

Figure  C-l. 
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4 

4 

-6C.8 

.000«-  *00 

.0001  »  1 

64.9 

•  0 

.0 

4 

6 

73?. 

*00 

.OOOF  ♦  <  ■’ 

169.8 

•  0 

.0 

4 

7 

87.8 

.OOOf *00 

•  OOOf  *o/» 

11.4 

.0 

.0 

4 

e 

-64.0 

.ooof*oo 

.OOOF  *00 

174.8 

•  0 

.0 

4 

9 

107. 

.000F ♦ 00 

.OOOF  *00 

174.  7 

•  0 

•  0 

2 

l 

-774. 

.OOOF  *00 

.OOOF  *00 

1  77.  7 

.0 

•  0 

2 

2 

-27.7 

,000F  *00 

.OOOF  *00 

88.0 

.0 

•  0 

2 

3 

201. 

.OOOF  *00 

.OOOF  *00 

36.3 

.0 

•  0 

2 

4 

33.7 

. OOOF  *00 

•  OOOF  *00 

116.8 

.0 

.0 

|  2 

4 

68.0 

.0008  *00 

.OOOF  *00 

144.6 

.0 

•  0 

2 

6 

-40.8 

. OOOF  *00 

•OOOF *00 

7.9 

.0 

.0 

2 

7 

-89.3 

.OOOF ♦ 00 

•  OOOF  *00 

94.7 

•  0 

•  0 

2 

6 

-18.0 

.OOOF  *00 

•  OOOF  *00 

6.4 

.0 

•  0 

2 

9 

18. 0 

.OOOF  *00 

•  OOOF  *00 

14.6 

•  0 

•  0 

3 

l 

979. 

.OOOF  *00 

. OOOf *00 

l  1.  1 

•  0 

.0 

3 

2 

-16.9 

.OOOF  *00 

•  OOOF  *00 

47.6 

.0 

•  0 

3 

3 

184. 

.OOOF  *00 

•  OOOF  *00 

42.8 

•  0 

.0 

3 

4 

-81.6 

.OOOF  *00 

•  OOOF  *00 

91.8 

.0 

.0 

3 

4 

64.4 

, OOOF  *00 

•  OOOF  *00 

43.9 

.0 

•  0 

3 

6 

-68.2 

.OOOf *00 

•  OOOF  *00 

44.9 

•  0 

.0 

3 

7 

119. 

.OOOF  *00 

•  OOOF  *00 

60.6 

•  0 

•  0 

3 

6 

-61  .4 

.OOOF  *00 

•  OOOF  *00 

42.7 

.0 

.0 

3 

9 

.OOOF  *00 

•  OOOF  *00 

4  7.6 

•  0 

.0 

1 

1 

1*6. 

. OOOF  *00 

•  OOOF  *00 

118.6 

•  0 

.0 

l 

2 

4.86 

.OOOF  *00 

•  OOOF  *00 

112.0 

•  0 

•  0 

1 

3 

-24.0 

• OOOF ♦ 00 

•000E*00 

118.9 

.0 

•  0 

& 

4 

47.8 

•  OOOF  *00 

•oooe*oo 

111.8 

.0 

.0 

1 

5 

-114. 

.OOOF  *00 

•  0  OOF  *00 

82.4 

•  0 

.0 

1 
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-?1.3 

.0006*00 

.0006 

♦00 

188.4 

.0 

.0 

7 

0 

41.8 

.0006*00 

.OOOE 

♦oo 

181.0 

.0 

.0 

7 

9 

68.4 

.0006*00 

.0006 

♦  00 

138.2 

.0 

.0 

5 

1 

246. 

.0006*00 

.OOOE 

*00 

31.0 

.0 

•  0 

5 

2 

-64.0 

.0006*00 

•  OOUf 

♦  00 

7.6 

.0 

.0 

5 

3 

-96.  0 

.0006*00 

.0006 

♦  CO 

6.  7 

•  0 

•  0 

4 

-103. 

.0006*00 

•  OOOF 

♦  00 

160.0 

.0 

.0 

5 

0 

-19.4 

.0006  *00 

•  OOOE 

♦  00 

147.8 

.0 

•  0 

5 

6 

144. 

.0006  *00 

.0006 

♦  00 

1  73.6 

.0 

.0 

* 

7 

-84.8 

.0006  *00 

.0006 

♦  00 

l  74.2 

.0 

.0 

5 

8 

-72.3 

.0006*00 

.0006 

*00 

.8 

.0 

•  0 

5 

9 

-96.7 

.0006*00 

.0006 

♦  00 

•  1 

.0 

•  0 

6 

l 

477. 

.0006  *00 

.000* 

*00 

81.7 

•  0 

.0 

Figure  C-l .  Program  Matrix  Listing  (Continued). 
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6 

7 

28,  3 

.OOOF  *00 

. OOOF  *  0 

1  76.9 

.0 

.0 

6 

3 

-33.8 

.0001 *00 

.00<*r 

106,  3 

.0 

.0 

6 

4 

-42.2 

•  OOOf *00 

.OOOF  *00 

109.  1 

.0 

.0 

6 

5 

4.72 

.oooe  *oo 

.OOOF  *00 

29.6 

.0 

.0 

6 

6 

-41.1 

•  000E  *00 

•  OOOF  *00 

29.  1 

.0 

.0 

6 

7 

-188. 

,0006  *00 

.OOOE  *00 

1  70.1 

.0 

.0 

6 

8 

-24.2 

, OOOF  *00 

.OOOF  *00 

42.6 

.0 

.0 

6 

9 

-76.0 

.OOOf  *00 

.OOOf *00 

72.2 

.0 

.0 

7 

1 

432. 

.onof  *oo 

.OOOF  *00 

6.6 

.0 

.0 

7 

7 

-274. 

.OOOF  *00 

•  OOOF  *00 

28,1 

.0 

.0 

7 

3 

-464. 

.OOOF  *00 

•  OOOF  *00 

.  7 

.0 

7 

4 

448. 

.000F*(10 

.OOOF *00 

14.4 

.0 

.0 

7 

8 

-14.4 

,00  E  *00 

.OOOf  *00 

14.0 

.0 

.0 

7 

6 

-130. 

•000E ♦ CO 

.OOOF  *U0 

15.3 

.0 

.0 

7 

7 

462. 

•  000  E *0 

.OOOE  *00 

27.3 

.0 

.0 

7 

8 

-223. 

.OOOf  *00 

.ooor  *00 

1. 1 

.0 

.0 

7 

9 

420. 

.OOOf  *00 

•  OOOF  *00 

1 4  7  .  3 

.0 

.0 

3 

1 

364. 

.OOOF  *00 

•OOOF ♦ 00 

1  12.0 

•  0 

•  0  j 

3 

2 

7.72 

.OOOF  *00 

.  OOOE  *00 

18.3 

.0 

•  0 

3 

3 

63.  ft 

.OOOF  *00 

, OOOE  *00 

94.3 

.0 

•  0 

3 

4 

-88.3 

.OOOF  *00 

.OOOf  *00 

162.1 

•  0 

.0 

3 

-3.97 

.OOOF  *00 

•  OOOF  *00 

35.5 

.0 

.0 

3 

6 

14.4 

•  OOOF  *00 

.OOOE  *00 

34.9 

.0 

.0 

3 

7 

-46.0 

•  OOOF  *00 

•  OOOE  *00 

21.4 

.0 

.0 

3 

8 

20.6 

• OOOE *00 

.OOOE  *00 

50.0 

.0 

.0 

3 

9 

76.1 

.OOOF  *00 

•  OOOE  *00 

73.2 

•  0 

.0 

3 

84 

2 

168. 284 

2. 8208817 

5 

3 

4 

90 

0 

180. 711 

.6866900 

5 

4 

5 

104 

0 

204.  393 

.7893097 

5 

5 

6 

no 

8 

220.972 

4.8140106 

5 

6 

7 

119 

l 

238. 822 

1.0497084 

5 

7 

5 

139 

0 

279.764 

.4410013 

5 

5 

6 

144 

l 

291. 175 

l  .0753403 

5 

6 

7 

140 

0 

300.294 

,5929406 

5 

7 

3 

142 

4 

*Z*.8I9 

2 . 0444567 

5 

3 

M 

£  NO  OF 

RECORD 

OR  A 

l-ENOt 

Figure  C-l. 


Program  Matrix  Listing  (Concluded). 
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666.076 


“2.619 

-  .  .V  2 

l ,  f 

1.6)? 

3*6 

-.01  2 

.  1 

.  » 

-6. 310 

-  .  C  t,  7 

.  6  7  6 

-••77 

2.562 

-.016 

—  •*91 

2.916 

793 

**./•  76 

,  2*6 

.292 

.607 

-.101 

.806 

-.025 

732 

-  .  C  39 

.61  ? 

-  .  u  *)  H 

.162 

.008 

-.099 

GU  7 

-.061 

.0  26 

.  120 

-.065 

52.106 

2.66  5 

22  .901 

-6. *  31 

2.665 

.173 

1.635 

- .  100 

22.903 

1.605 

12.536 

95? 

-6.631 

-.100 

-2 . 9S  7 

6 . 7  h  2 

.909 

•  09Q 

1.108 

.190 

-6.970 

-.263 

-3. 799 

-1.6*1 

10. 783 

.580 

6.862 

-1.821 

-2.668 

-.117 

- 1 . 1 C  7 

.  173 

6.35  1 

.265 

2.625 

-.20  7 

the  eigenvalues  ACE: 

The  COMPARISON  OF  EIGENVALUES: 

UNLOADED  PREDICTED  ACTUAL 

EXPERIMENTAL  MASS-LUAOEQ  MASS-LOADED 

0AT*  using  experimental 

-  _ -  SQUARE  MODAL-,  HAT  A 

MATRIX 


169.637 

m.2M 

206.396 
206*  6 A  3 
—  -2J8.039 
238.212 
279.993 
651.669 
686.68? 
I. 


159.686 

1  ?  6  •  8  6  6 

197.687 
206.699 

-233.10] _ 

2  38. 1  A* 

2  78.8  72 
650. 2B6 
683. 779 


139.976 
1 79. 832 
206.266 
206,685 

238.620 
238.660 
280.666 
6  S  3 . 66  8 
666.076 


.0393526 32761Q6 
-.175363006826 
.2073197098769 
37639?4,6??96  ] 


-•  16 70 1336*611  6 
.09  7  39a?  <,07  <,9  i  6 
.03701173771761 
“. 106329M6009S6 

•626l90O77S4«6 
- .065  706 86  766  369 


-.3585833139561 

1. 


r. 818857171 3678 
•02 66 *0 980*610 6 
. 161 69  12  m -0  3?  ) 
•112136026686? 


-.37666369  107*. 
•297093738?067 

-.QQ26C96026Q  ?2t6 
.06596301  37912  7 

l. 

•  15  3358768356  7 


|-.266G50  36  75  722 
j  —  •  1  10260  1  3  ?  1  0 9  3 
-.0006162736666787 


*■  i  i 

.  3  6t. 

0  *>9 

-  .  J6  7 

f-  1  l 

-.  16  ) 

7  6? 

.231 

66  2 

.00) 

2  1  i 

-1.097 

2  1  7 

-  .  J0  1 

u  vj 

1  16 

6  v9 

-.oul 

9  09 

-6.920 

0  9.) 

-.26) 

108 

- i. 789 

19.) 

“1.6*1 

At  i 

-.865 

H  66 

7.531 

0  lo 

-.173 

106 

.  6  6  5 

512 

-1.961 

.  '  1  J 

-.251 

-  .  -  ■«  1 

-.022 

-  .  J  t  .j 

-  .  060 

.1  2  7 

•  lU6 

-.115 

-.030 

-1  .  1  2s 

-.  185 

.13/ 

-.081 

-.  J91 

.  0?6 

,U69 

.  001 

lj. 73 j 

-2 . 668 

•  5  .0  u 

-.117 

6.8  6? 

-1.107 

-1.821 

.173 

J  10 

136 

-.  1  73 

.665 

2.69  7 

-.656 

-.  6  *>t» 

.  156 

.616 

-.270 

I  l  6 

•  G  to 
-.021 
-  •  2  Go 

.152 
.668 
.03  2 
-.015 
.  008 
6.363 
.266 
2.626 
-.  207 
.512 
-1.961 
.616 
-.270 

•  906 


PERCE  NT 

SQUARE 

PREDICTED 

PERCENT 

SQUARE 

error 

ROOT 

mass-loaoed 

Error 

ROOT 

predicted 

OF 

USING 

PR E D I  C  TE  0 

OF 

versus 

D 1  E  F  *  S 

PSEUDO  INVERSE 

VERSUS 

01 FF  •  S 

ACTUAL 

S  0*0 

1 9«0wS , 0COL S  > 

actual 

SO  *D 

-13.92 

.  39 

169.637 

-21,05 

2,  59 

1.82 

•  8  6 

18  1.281 

79 

2.65 

3.21 

.39 

206. 396 

-1.06 

1.66 

-.90 

.62 

206.663 

-.97 

1.69 

-  2.32 

3.87 

238.039 

.25 

3.80 

.22 

1.31 

238.212 

.  19 

3.27 

.63 

.06 

279,993 

.23 

1.26 

.75 

8.65 

661 . 669 

.  66 

8.71 

-3.80 

.72 

0.000 

100.00 

1.79 

Figure  C-2 


Edited  Output  From  Program  Matrix 


Sample  output 


»ri  o 


f  rom  LJISiSI.A  jr  «phics 


wOR«  a<M-r  |  s'  !  ns 


Mns<is*  l‘.  ,j(. 

Yinavjs- 

2  3D  ax i s  «  1 s • - o 

!N  A-S.  )-n  UNITS 


V! FWPOfNT 


XVU»-2.C0C£*ul 
y vu *  7*rocr*iw 
z  vu  •  S.00QM01 

IN  APS.  3-D  UNITS 


GRAPH  SFT-UP  I  GRAF30  » 


.ORIGIN 


X10nRlClN--2.000£*0l 
Y300» ICIN--1.00CI  *00 
Z3D0RIGIN—  S.OGuE-Ol 

STEP  SIZE  _ 


X  IDS  TP  •  1 • OuOE  +CI 
YIDSTP-  1 . OOOE ♦ 0  l 
Z3DSTP-  1 «CoOE «  0 1 


*3P*4X«  5.000£*01 

yid-ax-  s.oooe*oi 

l  3D’*  A  X  ■  S  .  OOQ  L  *0  1 


LOCATION  OF  CURRENT  PHYSICAL  ORIGIN 
X-  .SO  Y-  1.1?  INCHES 

FRO*  LOWER  LEFT  CORNER  OF  PAGE 


Figure  C-2.  Edited  Output  From  Program  Matrix 


(Continued ) 


L03 


^7h.  1*? 

•»'M .  *19 
323.2*5 
401  .*>7  8 
434.504 
456.524 

2  «  406 
.137 
1.409 
1.088 
-.005 
1.213 
-.287 
.  1C9 
-  -  .*32 

17.393 
-.84  3 
-1.041 
.45  3 
-.450 

-  - *6.72 Z 

2.2  16 
-.010 
1.229 


-.151 

.16  7 

1...W 

—  •  6  7  *, 

1 .4  74 

.  2  n  .1 

.  4  79 

6  l  0 

•  P  16 

-.071 

-.036 

-•  vJt>f> 

?  7H 

.^24 

-  •  u  3  3 

•  u2* 

.058 

-.442 

.  667 

-.00  7 

.00  3 

•  52  4 

•  o  12 

-.163 

-.01  * 

-.  3  4.) 

1.605 

•  1  5  \j 

.  i  *7 

.49} 

.  049 

-.160 

-.041 

.045 

-.116 

.2  02 

.002 

.  J  19 

-.032 

•  00  8 

-.128 

-.151 

-2.930 

.061 

-  2  .  J  1 9 

.‘-71 

-.217 

-.0  37 

.04  4 

-.162 

.081 

-.37i, 

-.1  39 

•  Ou  7 

.  1 18 

-.2  34 

.024 

-  .  5 «, 

.079 

-.  115 

164 

.005 

.004 

.001 

-.001 

—  •  1 6  9 

.667 

-.042 

.8  79 

.454 

•  lid 

-.  320 

-.84  1 

-1.041 

•  44  3 

-.450 

-6.227 

2 .736 

-.010 

1  •  229 

•  1  32 

*2  22 

-.013 

-.0  75 

.  1  75 

—  •  C  4  1 

.065 

-.098 

•  22  2 

l  .611 

l.  6  19 

-  .o67 

358 

-.4  73 

.073 

.962 

-.013 

1.639 

7.862 

.0  ?8 

-1.879 

-.349 

-.040 

1.594 

-.075 

-  .067 

.076 

.294 

.943 

-,  40l> 

-.132 

.022 

.175 

- 36. i 

*1.879 

.941 

9.367 

-1.835 

-.307 

-1.267 

-.041 

-.473 

349 

-.400 

-1.835 

1.257 

.  159 

-.047 

•  C  6  6 

.073 

—  •  04  U 

-.132 

-.307 

.  159 

•  072 

-.048 

-.098 

.962 

1 . 594 

.022 

-1.247 

-.047 

-.048 

1.214 

5  ABF  : 

—  UNLOADED  —  -PfiECICUQ 
EXPFPr"FNTAl  8ASS-LPA0FD 
DAT4  USING 

5  OU  APE  800AL 
MATRIX 

137.275  137.085 


137.275 

- 207.085 

239.182 
280.55  3 
294.848 
327.081 
394.224 
- *14.164- 

453.685 


ACTUAL 
"ASS-LOAOf 0 
EXPERI CENTAL 
OATA 

149.8  72 


-  18  X*  200 - 202.59a 

227.943  238.415 


260. 786 
286.278 
327.150 
393.434 
— 412.494 
453.491 


276.142 
301.519 
32  3.245 
401.578 

- 434.504 

456.524 


.008457103552768 
-.1817037949291 
-.1973675842933 
-.04608525837128 
.22681980  3*452 
.1959742317667 
.01445647107925 
-.1195120697267 
.0389728  2641439 
-.0513527429707 
.2057643600871 
1  . 

.5553119430^81 
-.  7256695  74  712  0 
-.129154047651 1 
-.07121916674848 
.2303  7 7e  71  7*47 
-.32416158678)1 
•  092  65044P 7  >9^ 
-.200521216  74  !  1 
.7396044  7  77  76  7 
-.441  37 60 8^*477 


PERCENT 

SQUAR  £ 

PREDICTED 

PERCENT 

S  QUA  ft | 

error 

ROOT 

^ AS  S-LOaDEO 

ERROR 

boo  r 

’REDICTFO 

OF 

USING 

PREOICTEO 

OF 

VERSUS 

0 1  F  F  •  S 

PSEUDO  INVERSE 

VFRSUS 

D1FF  •« 

ACTUAL 

SO  *0 

19R0*S%8CQLS» 

actual 

SQ  1  D 

8.53 

.46 

137.275 

8.41 

.43 

10.56 

-  *44 

-  -  207*085  ..  . 

-  -2.21 

2.67 

4.39 

3.96 

239.  182 

32 

4.28 

5.56 

1.85 

280. 553 

-1.60 

3.  04 

5.05 

1.62 

294.848 

2.21 

2.53 

-1.21 

.77 

327.881 

-1.43 

•  89 

2.03 

.17 

394.224 

1.83 

.  35 

4*74  - 

- *  16  - 

—  414. L64 

4.68 

.26 

.66 

4.08 

.001 

100.00 

2.99 

Figure  02.  Edited  Output  From  Program  Matrix  (Continued) 
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2  79.  764 
2*3. 1 76 
300.294 
324.819 


-.678 

7  .  “66 

.  ?  3  7 

l  .  V  6  s 

.  >  4C 

—  .SCO 

1.9*7 

.  fS  . 

-  .  >9  9 

•  9  4  l 

66#. 

-.749 

1.814 

•  77  7 

-.117 

-.277 

-.6  96 

-.111 

.  6  n* 

.71* 

1.917 

.  1  7  6 

-.  7  7« 

•  0*3 

.2  36 

1  4  •> 

-.961 

.0  »8 

-.107 

1.16  7 

.071 

-  .  P  >  M 

-.  1  IN 

1  •  9i)M 

.113 

1  .09  4 

—  •  L'  l  9 

.  1  79 

•  C  80 

-.828 

-.  197 

779 

.  isi 

t  7/: 

.1  87 

-  .  M  *} 

-.611 

.362 

1.610 

1.162 

.2  16 

-1.229 

7. *01 

i  io 

-1  .  i  M 

-t  .  9 NO 

.08  1 

-1.741 

-.13C 

.702 

1.0  72 

-.62  9 

-.Oil 

-.887 

-1.171 

!  .072 

2.19  7 

-.439 

-.046 

-1.093 

-1.950 

-.629 

4  19 

2.986 

-.066 

.18? 

.081 

-.OJ  l 

—  .U46 

-.066 

.0  36 

•  004 

-  t  •  2  4  1 

-.88? 

-1.091 

.  187 

•  004 

2.561 

.877 

-.292 

-1 .169 

,08b 

.1*  l  7 

—  •  6  9  ? 

-.064 

.645 

.979 

-.557 

.005 

-.921 

.*'9  .19?  .0?.’ 

EICENVJIUFS  ARf: 

CP"f’AS(SnM  rtf  FlCFHVllU FS: 

l  .010 

-.069 

-1.616 

.10  0 
.  326 
.  -i  *>5 
.7b  l 
.022 
.  7  07 
.  714 
.765 

.  4 22 
.877 
.292 

.  169 
.0  86 
.017 
.•392 

.803 

.184 

.9S4 


.96" 
.001 
-.021 
.12  8 
-.064 
-.  764 
•  l  78 

.009 

.74  7 
-.064 
.645 
.  9  79 
-.557 
.006 
721 

184 

.621 

.241 


1.04  1 
.  166 

.  261 
.041 
-.(  *9 
-1.6  >0 
•  2  M 
.181 
1.544 
.458 
.  1  42 
.072 
1.010 
-  *  0  6  9 
-1.535 
.954 
.241 

1.849 


unloaded 
EXPEB I Mf NT AL 
DATA 

167.266 

180.424 

207.2*8 

228.543 

238.610 

279.684 

299.232 

303.621 

328.621 

•07756036581111 
.569376  f>C96?6? 
~.41941386R621«- 

•  3240401 302842 

1. 

-.3249416654984 

•1551939140735 

-.1024288017247 

•03645737474771 

•3779360199776 

1. 

-.5596484697106 

•  464  7784  I  794 1 
.01854646*0  1*«, 

34*  1649996014 

•  197241 711794 1 
14784  7 6367  76 

.2729614168009 

•  56  73664094607 

-.830968074780#, 

. 787559900419 
1147410716766 
.1441769077 77? 


PREDICTED 

"45  S-Ln aof  r» 

USING 

S0U4P  E  *OPM. 
"4 T a  I  x 
105.091 
179.65b 
199.139 
227.743 
237.447 

262.265 
290.583 
301.309 
125. 109 


Af  TUAL 

pf »CI NT 

SOU  A®  f 

PREDICTED 

PERCENT 

SQUARE 

ass-loach  o 

f  R ana 

ppn  r 

*AS  S-LUADFO 

f  BROB 

ROOT 

XPf  R  r  **£  N  T  AL 

PREDICTED 

Of 

USING 

PREOIClf D 

cc 

04  T  A 

VERSUS 

01  FT  *S 

PSfUOOjNVEBSf 

Vf  asus 

01 ff  *s 

ACTUAL 

SO  *  D 

I9B0WS.8C0LS1 

ACTUAL 

SO  #D 

1  68.7  85 

37.55 

1.47 

167.766 

.61 

1.63 

180.711 

.58 

2.14 

180.424 

.16 

2.21 

2  u 8 . 19  3 

4.34 

1.26 

20  7.288 

.63 

2.  16 

220.972 

-3.06 

i.  61 

228.643 

-3.41 

1.60 

7  38.822 

.68 

1  .64 

218.610 

•  12 

2.17 

2  79,  764 

6.76 

1.4? 

279.664 

.03 

1.91 

793.1 75 

.88 

7.49 

299.232 

-2.07 

1.56 

K 0.294 

-.34 

7.62 

303.621 

-l.il 

2.66 

174.819 

-.16 

1.68 

.002 

100.00 

1.60 

Figure  C-2.  Edited  Output  From  Program  Matrix  (Concluded) 
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Figure  C-2.  Sample  NAS  TRAN  Deck  (Continued)  . 
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Figure  C-2.  Sample  NASTRAN  Deck  (Concluded). 
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Fiqure  C-3.  Data  Acquisition  Figure  C-3.  (Concluded). 
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Frank  Broderick  Atkinson  was  born  on  13  November 
1947  in  Nashville,  Tennessee  to  Robert  T.  and  Ruth  B. 
Atkinson.  After  graduation  from  Irving  Senior  High 
School  in  1965  he  attended  Arlington  State  College 
(now  the  University  of  Texas  at  Arlington) .  During 
his  studies  there  he  participated  in  the  cooperative 
education  program  as  a  Weight  Control  Engineer  at  Bell 
Helicopter  Company  in  Hurst,  Texas.  Following  gradua¬ 
tion  in  1970  with  a  Bachelor  of  Science  degree  in  Aero¬ 
space  Engineering,  he  enlisted  in  the  U.  S.  Air  Force  as 
a  Ground  Radio  Communications  Equipment  Repairman.  He 
subsequently  attained  the  rank  of  Sergeant  whereupon  he 
was  selected  to  attend  Officer's  Training  School  (OTS) 
at  Lackland  AFB,  Texas.  Upon  graduation  from  OTS  in 
January,  1975,  he  was  assigned  to  the  3246  Test  Wing, 

Guns  and  Fuzes  Division,  Eglin  AFB,  Florida  as  a  Wing 
Munitions  Test  Engineer.  During  his  tenure  at  Eqlin 
AFB,  he  was  selected  three  times  as  his  division’s 
nominee  for  the  Directorate  of  Test  Engineering  Test 
Engineer  of  the  Quarter  Award,  He  also  represented  the 
Directorate  of  Test  Engineering  as  the  nominee  for  the 
Lt.  Robert  L.  Sullivan  Award  for  the  most  outstanding 
junior  officer.  Capt.  Atkinson  was  assigned  to  the 
School  of  Engineering  of  the  Air  Force  Institute  of 
Technology  in  June  1979  in  the  Graduate  Aeronautical  Engi¬ 
neering  Program. 
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